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We review recent theoretical developments in the study of the structure of jets 
that are produced in ultra relativistic heavy ion collisions. The core of the review 
focusses on the dynamics of the parton cascade that is induced by the interactions 
of a fast parton crossing a quark-gluon plasma. We recall the basic mechanisms 
responsible for medium induced radiation, underline the rapid disappearance of 
coherence effects, and the ensuing probabilistic nature of the medium induced 
cascade. We discuss how large radiative corrections modify the classical picture 
of the gluon cascade, and how these can be absorbed in a renormalization of the 
jet quenching parameter q. Then, we analyze the (wave)-turbulent transport of 
energy along the medium induced cascade, and point out the main characteristics 
of the angular structure of such a cascade. Finally, color decoherence of the in¬ 
cone jet structure is discussed. Modest contact with phenomenology is presented 
towards the end of the review. 
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1. Introduction 

A jet is a cone-shaped beam of energetic particles which originate from a hard 
parton produced in a high energy collision. Jets have played an important role in the 
development of Qauntum Chromodynamics (QCD), with three jet events in e + e - 
annihilation providing a direct experimental evidence for the existence of the gluon, 
and the quark-gluon vertex. 1,2 Moreover, perturbative QCD correctly describes 
the soft hadronic activity surrounding jets. For instance, the inclusive intrajet 
hadron distribution, the so-called fragmentation function, exhibits a humpbacked 
shape that reflects the suppression of soft gluon radiation at large angle as a result 
of destructive interferences, as predicted by QCD. 3 “ 5 Today, jet technology has 
become a precision tool embodied in powerful Monte Carlo event generators, 6 ^ 8 
and involve elaborate perturbative QCD calculations (culminating at the NNLO 
level), that are not only relevant for further testing strong and weak interactions 
at high energy, but also in searching for physics beyond the Standard Model at the 
LHC (see e.g. Refs. 9,10 recent reviews). 

In the context of heavy ion collisions, jet occupy a very special place. They are 
produced in hard collisions over tiny space-time regions, and the partons that orig¬ 
inate form such hard collisions propagate through the hot matter produced in the 
bulk of the heavy ion collisions. As such they constitute ideal “test particles” whose 
study could reveal the properties of the medium they traverse before hadronizing. It 
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was recognized long ago, 11 for instance, that quark and gluon jets may lose energy 
while crossing a quark-gluon plasma produced in heavy ion collisions, resulting in a 
strong suppression of high-p t hadrons. This phenomenon, commonly referred to as 
jet-quenching , was in fact observed at RHIC 12 15 and at the LHC, 16-18 providing di¬ 
rect evidence for the creation of dense QCD mater in heavy ion collisions. The high 
energy of the LHC make jets a very powerful diagnostic tool: they are produced in 
large number, have large energies (with p t ranging from 100 to 300 GeV, an order of 
magnitude larger than at RHIC). Furthermore, the detector capabilities of ALICE, 
ATLAS and CMS, allow for the exploration of a variety of new and more detailed 
observables, among which, the jet nuclear modification factor, the dijet asymmetry, 
the jet profiles. 16, 18-27 All this is opening new perspectives in the field and invites 
new theoretical efforts. 

The theory of jet quenching in heavy ion collisions started developing in the 
early 90’s with the first studies of medium induced gluon bremsstrahlung by Gyu- 
lassy and Wang, 28,29 accompanied by those of Baier, Dokshitzer, Peign, Mueller, 
Schiff 30-32 and Zakharov, 33,34 where it was realized that radiative parton energy 
loss is the dominant mechanism for jet quenching in large QCD media. This pi¬ 
oneering works were followed by many others emphasizing various aspects of the 
phenomena and using different formalisms. For instance, assuming the medium 
to be dilute, Gyulassy, Levai and Vitev 35,36 (see also Ref. 37 ) developed a formal¬ 
ism in which the interactions with the plasma can be calculated order by order in 
opacity. In a similar approach, due to Guo and Wang, the so-called High-Twist 
(HT), 38,39 the first medium correction to the the leading order splitting functions 
is computed and embedded in a medium-modified version of the DGLAP equation. 
Considering the plasma to be dense, an approximation of multiple soft collisions 
was applied by Armesto, Salgado and Wiedemann 40,41 where multiple gluon radi¬ 
ations are treated as independent and described with a Poisson distribution. 42 An 
analogous approach, accounting for thermal effects, is due to Arnold, Moore and 
Yaffe. 43-45 

While based on similar assumptions (see Refs. 46-48 for reviews), the various ap¬ 
proaches fail to provide a consensual quantitative description of the RHIC data. 49 
There are perhaps two reasons for that. First, the relatively low momentum range 
accessible at RHIC, of the order of 10 GeV, may question the applicability of weak 
coupling techniques. Second, each of the models has limitations related to ap¬ 
proximations done regarding the nature of the medium and the kinematics. For 
instance, none of the models addresses from first principles the issue of in-medium 
multiple parton branchings and resulting interferences, and thus rely on strong as¬ 
sumptions on the space-time structure of the in-medium parton shower. Moreover, 
at RHIC, scenarii based on radiative energy loss proved to be inconsistent with the 
strong suppression of the single non-photonic electrons resulting from semi-leptonic 
decays of hadrons carrying heavy quarks, 50 reviving the interest for collisional en¬ 
ergy loss. 51-54 At the same time, evidence for almost perfect fluid behavior of the 
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quark-gluon plasma, 55-5 ' diverted the attention towards strong coupling techniques, 
providing an alternative perspective on jet quenching based on the AdS/CFT cor¬ 
respondence (see Ref. 58 for a recent review, and references therein). 

In parallel to the aforementioned developments, the theory of jet-quenching 
based on essentially weak techniques of QCD, has progressed within the last few 
years. The set goals are the understanding of the structure of jets in matter, and 
how this structure is modified with respect to that of well understood jets in vac¬ 
uum 59-72 (see Ref.' 3 for a recent review), and the calculation of higher-order correc¬ 
tions, 74-80 addressing some of the limitations of the standard approaches of energy 
loss. Among the major results of these efforts, it was shown that in the limit of 
asymptotically large or dense media successive medium-induce branchings are quasi¬ 
local and independent from one another, 65,78 providing a theoretical basis for the 
approximations made in some models of jet quenching. 40,42,45 Color coherence, that 
characterizes jet structure in vacuum, has been shown to be significantly altered in 
a QCD medium. 59-61 As a result, medium-induced color decoherence suppresses 
the interferences between successive gluon radiations. There is also an effort to 
systematize the calculations of multiple parton branching in the framework of Soft- 
Collinear Effective Theory (SCET). 76,81 Last but not least, radiative corrections to 
some jet-quenching observables were recently carried out, 74,75,7 ' - ' 9,82 in addition to 
calculations of the jet quenching parameter on the lattice. 83,84 The present review 
will present a summary of some of these developments. 

We should also mention that the theoretical developments mentioned above are 
accompanied by ongoing efforts to develop Monte-Carlo event generators, that en¬ 
compass some of the recent advances in the microscopic physics, as well as a realistic 
modeling of the geometry of the plasma and its hydrodynamical evolution, in or¬ 
der to carry out systematic studies of jet observables with all sorts of experimental 
biases, such as, p t cuts, triggers, etc. The most notable ones are: JEWEL, 85,86 
Q- PYTHIA 8 ' and MARTINI, 88,89 that account for the BDMPS-Z radiation mech¬ 
anism, in addition to PYQUEN/HYDJET, 90 YAJEM 91 and MATTER. 92 These 
developments will not be covered in the review. Also, the phenomenology of jet 
the quenching is beyond the scope of this review and will be covered in another 
contribution to the present volume. 93 

Before we embark on the theoretical discussion, a word of warning. For the sake 
of simplicity we shall most of the time deal with the ideal scenario in which the 
parton cascade is purely gluonic and interacts with a uniform and static plasma. 
We also assume the medium to be dense enough so that we can apply the multiple 
soft scattering approximation. Corrections to this idealized picture can be calcu¬ 
lated, and this will be mentioned at proper places, but those corrections do not 
change the qualitative picture. The review is organized as follows. In Section 2, the 
physics of parton branching in vacuum and the BDMPS-Z mechanism for medium- 
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induced gluon radiation are introduced. In Section 3, we generalize the BDMPS-Z 
spectrum to multiple branchings and develop a probabilistic picture of in-medium 
gluon cascades. Then, in Section 4, we analyze radiative corrections to momen¬ 
tum broadening and parton energy loss and show that they can be reabsorbed in 
a renormalization of the quenching parameter. A complementary discussion on 
NLO corrections to energy loss and p ±-broadening in high temperature plasma is 
presented in this volume Ref. 94 : note however that these corrections are related 
to soft plasma modes and are of a different nature than the radiative corrections 
involved in the renormalization of q. In Section 5, we examine the time evolution 
of the angular and energy distribution along the medium-induced cascade, away 
from the jet core. In Section 6, we discuss the physics of in-medium jet decoherence 
that affects the inner jet structure. Finally, Section 7 presents a brief summary and 
outlook. 


2. Jets and parton cascades of two kinds 

Perturbative QCD successfully predicts the emergence of jets in the final state of 
high energy collisions such as e + e _ annihilation and hadronic collisions. As a matter 
of fact, the collimation of jet particles relates directly to the collinear nature of QCD 
parton splittings. The underlying mechanism is akin to QED bremsstrahlung: an 
electric charge that suffers an acceleration subsequently radiates photons collinearly 
to its direction of motion. Likewise, quarks (and gluons) carry a color charge and 
thus can act as sources of collinear gluon bremsstrahlung. Besides the confining 
forces of strong interactions that eventually turn patrons into observable hadrons, 
the difference with QED (at the partonic level) lies in the fact that, contrary to 
photons, gluons may also radiate affecting the general pattern of the collimated 
cascade. The nascent parton possesses typically a large virtuality Q 2 , and therefore 
radiates gluon, which themselves radiate other gluons, forming a so-called parton 
cascade. Because this cascade generically propagates in vacuum, we shall refer to it 
as a vacuum cascade. Such a cascade is accompanied by a degradation of the virtu¬ 
ality Q, which can be used to define an evolution “time”, t = 1 yiQ 2 /Qq, where Q o is 
the scale that marks the end point of the parton cascade where hadronization takes 
place. Typically, Q o ~ Aqcd- Hadronization is a non perturbative phenomenon 
and one may fear that it erases all memory of the detailed structure of the preced¬ 
ing partonic phase. However, this is not so. It turns out that inclusive observables 
are to a large extent determined by the parton dynamics prior to hadronization. 
This remarkable property, known as the local parton-hadron duality, 3,4 is largely 
responsible for the success of perturbative QCD in describing jet observables to a 
high degree of precision. 

When the energetic parton propagates in a medium, another mechanism may 
generate radiation, beyond that due to the possible virtuality of the parton: dur¬ 
ing its propagation, the parton collides with the constituents of the medium, and 
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each of these collisions may in principle trigger radiation. In fact, because of this 
mechanism, even an on-shell parton entering a medium will be able to radiate and 
generate a cascade. We shall refer to such a cascade as to a medium-induced cascade. 
We therefore encounter two types of parton cascades: 

(1) The vacuum cascade that forms after a hard process that creates a parton with 
a large virtuality. 

(2) The medium-induced parton cascade that forms due to multiple interactions of 
the jet with a QCD medium. 

This review will be mostly devoted to an analysis of the medium induced cascade. 
But in our discussion we shall often refer to the vacuum cascade in order to empha¬ 
size some of the distinctive properties of the medium induced cascade. 

In both cascades, coherence and interferences play an important role by sup¬ 
pressing radiation in particular regions of the spectrum. For the vacuum cascade, 
this is the case of gluons with large transverse wavelengths which do not resolve 
the individual color charges carried by the jet, and thus are radiated coherently off 
the total color charge of the jet. In the case of medium induced radiation, we have 
gluons with long formation times that do not “see” the individual constituents on 
which they scatter, and appear to be radiated coherently off many scattering cen¬ 
ters. In both cases the gluon bremsstrahlung spectrum is suppressed with respect 
to the incoherent case. We now discuss these phenomena in some more detail. 

2.1. Jet structure in vacuum: transverse color coherence 

Experimentally, a jet is defined by a collection of particles that are contained in a 
cone of given opening angle Oj e t- The characteristic hard scale for the perturbative 
evolution is then Q = E0j et Aqcd- It corresponds to the largest transverse scale 
(with respect to the jet axis) that can be carried by one of the jet constituents and 
serves as an upper bound for the jet invariant transverse mass. The phase space 
between the hard scale Q and the non-perturbative scale Q o ~ Aqcd is filled by 
the vacuum cascade. 

The microscopic mechanism generating the cascade is bremssthalung. It is com¬ 
mon to both the vacuum and in-medium shower. In the collinear limit, the branch¬ 
ing probability into two offsprings, with energy fractions z and \ — z respectively, 
reads 

p{z)Az ' {2A) 

where a s is the strong coupling constant, and P(z) the Altarelli-Parisi splitting 
function. 95 The expression (2.1) of the branching probability exhibits a collinear 
divergence that corresponds to a mass singularity, and a soft divergence, when 
k± —> 0 and z —> 0 (z —>• 1), respectively. The latter is cured in inclusive quanti¬ 
ties by real-virtual cancellations, while the former is regularized by cutting off the 
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transverse momentum integration at the scale Qq. The collinear divergence pro¬ 
vides a large logarithm IuQ 2 /Qq that may compensate the small coupling constant 
in the expression (2.1) of the emission probability. When alnQ 2 /Q\ ~ 1, multiple 
branchings need to be resummed. In this resummation, successive branchings are 
strongly ordered in their transverse momenta, with the largest momentum emitted 
first. 9 5-9 7 

It was soon realized that this picture is not complete and fails to describe the 
soft sector of jet observables. This is because it does not fully take into account the 
effects of color coherence: The radiation pattern of a collection of color charges that 
form a jet, as they originate from a single parton, does not simply reduce to the 
incoherent sum of the radiation spectra of the individual charges. If the transverse 
wave length A_l ~ l/fcj_ of the radiated gluon is larger than the transverse size of 
the jet rj_ ~ 0j e tt/, at the time tf = uj/k\ when the gluon forms, the latter does 
not resolve the fine structure of the jet. The radiation is then triggered instead by 
the total charge of the jet, i.e., the color charge of the original parton. This leads to 
a suppression of the radiation as compared to that expected form the various color 
charges contained in the jet. 

Such considerations, including the fact that gluons carry a color charge and 
thus can radiate as well, 98-101 and accounting also for for energy conservation, 
constitute the basis of the jet algorithms in vacuum that were first introduced in 
the Marchesini-Webber model. 102 Along a QCD shower a strict ordering in angles 
accounts for the color coherence and the evolution equation for the inclusive parton 
distribution reads, 


dD(x, Q) 
dHQ 2 /Ql) 



dz a s {k±) 
z 27r 


P{z)D{x/z, zQ). 


( 2 . 2 ) 


where k± = z( 1 — z)Q and P(z) stands for the regularized Altarelli-Parisi splitting 
function 95 (the dependence on different parton species is omitted). Eq. (2.2) is 
referred to as the Modified-Leading-Log-Approximation (MLLA) 3,100 " 11 ^ 1 as it differs 
from DGLAP 95-9 ' by the fact that in the r.h.s. the parton distribution scales with 
zQ rather than Q to account for angular ordering. Eq. (2.2) describes in reality 
the evolution of the parton shower as a function of the jet-opening angle @j e t with 
Q = 0jet-E and E the energy of the parent parton: The infinitesimal variation of the 
parton distribution is given by an initial splitting of the parent parton generating 
a subjet characterized by the energy E' = zE and the corresponding scale is Q' = 
E'Q jet = zQ. 


2.2. Medium-induced gluon radiation: longitudinal coherence 

In a medium of size L, an energetic parton experiences multiple collisions with the 
medium constituents. In principle each of these collisions puts the gluon off-shell, 
so that it may subsequently emit radiation, after a typical time, commonly referred 
to as the formation time , tf ~ u>/k^_, where ui and k± are the energy and the 
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transverse momentum of the radiated gluon. However, this simple picture may be 
strongly affected by coherence effects in multiple scattering. 

Such effects were studied long ago in the context of QED, by Landau, Pomer- 
anchuk and Migdal (LPM), 103,104 and are generically referred to as the Landau- 
Pomeranchuk-Migdal (LPM) effect. The QCD analog of the LPM effect 0, was inves¬ 
tigated by Gyulassy and Wang (GW) 28,29 first, then by Baier, Dokshitzer, Mueller, 
Peign, Schiff and Zakharov (BDMPS-Z). 30-34 The main focus of these early studies 
was the energy loss caused by induced gluon radiation. 42 For large media, radiative 
energy loss scales as L 2 where L is the length of the matter traversed, as opposed 
to collisional energy loss that scales as L. 105 

We shall return extensively in this review to the BDMPS-Z mechanism, but 
it is useful at this stage to have an heuristic understanding of the basic physics 
involved. Consider then a hard quark or gluon of energy E propagating in a dense 
QCD medium. The collisions of the parton with the medium constituents induce 
radiation, with the typical formation time of a radiated gluon given by tf ~ w/fc 2 . 
Now, during that time scale, the gluon fluctuation experiences multiple kicks, which 
may cause its transverse momentum to increase by an amount k 2 ^_ ~ qtf where q is 
a diffusion coefficient to be specified shortly. Assuming that these multiple kicks are 
the only source of transverse momentum, and combining the two relations above 
one finds that the typical time for induced emission, to which we shall refer here as 
the branching time is given by 



(2.3) 


Multiple collisions cause the broadening of the transverse momentum. In the 
regime where the collisions are soft, and where a large number of collisions are 
needed to change the transverse momentum by a significant amount, a diffusion ap¬ 
proximation becomes appropriate. In this case, it is customary to define a transport 
coefficient q that describes the diffusion in momentum space. Assuming that the 
dominant color field fluctuations of the medium are characterized by a scale inn, 
to which we refer as the Debye screening mass, the typical transverse momentum 
exchanged with the medium in a collision is ~ mo- Then the transverse momentum 
squared acquired by the hard parton traveling a distance L is, by definition of q, 
qL. Equivalently, if A e i denotes the elastic mean free path, we have m 2 D ~ qX c \. 
The coefficient q is referred to as the jet quenching parameter. 

Now, when Tb r < A e i, each of the N = L/X e \ collisions act incoherently and 
produce a maximal radiation intensity oc N. This is the so-called Bethe-Heitler 
regime which yields a spectrum 



d I _ L 


(2.4) 



On the other hand, when r^ r A e i, a number iV co h = 7br/A e i of the collisions 
a See Ref. 106 for a recent review on the LPM effect in QED and QCD. 










March 23, 2015 0:17 


World Scientific Review Volume - 9.75in x 6.5in 


QGP5-BM-FINAL 


J. -P. Blaizot and Y. Mehtar-Tani 


act coherently as a single collision to radiate one gluon. As a result the radiation 
spectrum scales as N/N coh and reads. 31-34 bc 

= ^ < 2 ' 5 ) 

where lo c is such that Tb r ( w c ) = L. When lo = lo c , ujdl/duj ~ a, as if we had 
a single collision. The full Bethe-Heitler result Lodl/dio ~ aN is recovered when 
lo = q \ 2 , = lo B h- 

Note that for lo > lo c the spectrum is strongly suppressed. In this case, the 
radiated gluon cannot distinguish the hard event that produced its radiator from 
the multiple scattering of the latter with the medium. As a result, the suppression 
1 /y/uj turns into a stronger one, that is 1/w 2 . A complete calculation, accounting 
for finite size effects, with the two limiting cases above, yields 31 

lo^- ~ 2d In |cos (fiL)|, (2.6) 

where d = ol s Cr/tt, with Cr = Cf,a is the color charge of the source of radiation, 
and 


n = 



1 

T br (w)’ 


(2.7) 


By radiating gluons, an energetic parton, with energy E w c , that traverses a 
dense QCD medium, loses a mean energy, d 

r°° (\ t 7T 

A E = / d loco — = —dqL 2 , (2.8) 

Jo duj 4 


Indeed, since the spectrum Eq. (2.6) falls as 1 /yjw when w<w c (cf. Eq. (2.5)), 
the lo integral in (2.8) is dominated by gluon frequencies of the order of ui c . This 
formula relates the mean-energy loss to the typical transverse momentum squared 
accumulated by the gluons during the maximal formation time tf ~ L, (k±) ~ qL , 
that is, A E ~ (k\)L. 

Consider now the number of radiated gluons with frequency larger than lo, 

r u ° d I 

N(lo) = / da/ —. (2.9) 

Jui 


When lo <C lo c , the number of radiated gluons is dominated by the lower bound, lo, 
of the integral in Eq. (2.9), for the spectrum rises as 1/w 3 ^ 2 . Thus, we obtain 


N(lo) ~ 7 ^®. (2.10) 

V LO 

h I stands for “Intensity”, a dimensionless quantity proportional to the number of gluons. 
c In QED, the spectrum scales as OL S y/ajuj c /E. It depends explicitly on the energy of the projectile. 
At asymptotically high energies, the photon radiation dies out, while according to Eq. (2.5) gluon 
radiation does not. 

d For parton energies E <C oj c we find that A E ~ yJ~E 
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where 


uj s (L) = a 2 lu c (L)- 


( 2 - 11 ) 


As long as N(uj) is smaller than one it can be identified with a radiation probability, 
i.e., N(co) ~ P(w) < 1. When the average number of emitted gluons exceeds unity, 
higher orders in perturbation theory, that involve multiple gluon radiations, become 
relevant processes that need to be dealt with. According to Eq. (2.10), this regime 
is achieved when w < lo s (L). This criterion can be looked at from a different 
perspective. As we shall see later, the branching rate for a gluon with energy uj is 
given by 

= rrv (2 - 12) 

T br (w) f*(w) 

where, accordingly (see Eq. (2.3)) 





(2.13) 


is the time it takes a gluon with energy to to branch, a quantity that one could 
also refer to as an “inelastic mean-free path”. Then, the condition oj < u) s {L) is 
equivalent to the condition f*(w) < L , that is, during the time f*(w) a gluon with 
energy uj will branch at least once during the time it spends in the matter. When 
applied to the leading particle, the quantity f*(.E) will be referred to as the stopping 
time in Section 5. 


3. Probabilistic picture of medium-induced gluon cascade 

We start now examining more closely the properties of the medium-induced cascade 
and begin by putting in place various elements of the formalism that will be used 
throughout this review. 

3.1. Non relativistic dynamics in the transverse plane 

A high energy jet originates from an ultrarelativistic parton, whose direction of 
propagation is chosen as the z axis. We choose light-cone variables® such that the 
4-momentum p of a parton reads 

P + = \{E + p z ), p~=E-p z , and p = {p x ,p y ), (3.14) 

where E is the energy, p z the longitudinal momentum, and p the transverse mo¬ 
mentum. The angle 9 formed by the 3-momentum p and the jet axis is given by 
p_ l = \p\ sind. Considering that jets are collimated objects, and that their angular 
broadening is small, we assume in all our calculations that 9 <C 1, and accordingly, 

p + ~ E, and \p\ = p± ~ E6 , (3.15) 

e With this choice, the 4-coordinate vector is x = (x + ,x~ ,x), where = (t + z)/ 2, x~ = t — z, 
and x = ( x, y). 











March 23, 2015 0:17 


World Scientific Review Volume - 9.75in x 6.5in 


QGP5-BM-FINAL 


10 


J. -P. Blaizot and Y. Mehtar-Tani 


while the light-cone time can be identified with the real time, 

x + ~ t ~ (3.16) 

We work in the light-cone gauge A + = 0, a convenient choice in the present context. 
This is because the current J _ , which describes the target through which the jet 
propagates, is coupled to the -I- component of the gauge field, which vanishes in 
this gauge. In turn, gluon radiation off particles propagating in the — z direction is 
suppressed, and only the jet constituents, propagating in the +z direction, radiate. 
Furthermore, the light-cone gauge allows for a partonic picture of the radiation 
dynamics. 

The hard gluons that compose a collimated jet travel along nearly straight-line 
trajectories, and interact with the color field of the medium. Only the — component 
of this field couples to the high energy gluon with the helicity conserving vertex, 
V l; i + A~ ~ p + A~5 l G where i,j = 1,2 depict the two transverse components of the 
physical polarization vector e\. The transverse component of the field can be re¬ 
moved by making use of the additional degree of freedom of the light-cone gauge. 107 
Moreover, the gluon with very large longitudinal momentum p + , propagates near 
the light-cone branch, x~ = 0, so that, finally, the only component of the medium 
field that matters is A(t,x) = A~{x) | x - =0 - 

The choice of the light-cone gauge along with the eikonal approximation for the 
vertex imply that the gluon propagator in the medium background field is diagonal 
in helicity space, i.e., G 1 - 7 ~ <5 y G, with G obeying the following equation, 

[□ — igd + A~ (x) ■ T] G( x, x') = 5{x — x r ), (3.17) 


where g is the gauge coupling constant. We have already mentioned that the jet 
probes the medium gauge field only near the light cone, i.e. x~ ~ 0. We may 
then assume that over such small distance, the medium gauge field is constant, and 
perform a Fourier transform with respect to a; - . 1 We define 




i 

2E 


G ( X -a/), 

Z7T 


(3.18) 


and from Eq. (3.17) above, one finds that obeys the following 2-dimensional 

Schrodinger equation 


i 


d_ 

Ft 


di 

2 E 


+ gA(t,x) 


(x\g(t - t')\x’) 


iS(t — t')S(x — x'), 


(3.19) 


where the background field, A = A a T a , with T a the generators of SU(2 >) in the 
adjoint representation, is a matrix in color space, and so is the propagator Q . In 
the high energy limit, the interesting dynamics is therefore localized in the trans¬ 
verse plane. It is the 2-D equivalent of a non-relativistic particle in a background 

f Since x~ is the conjugate variable to ~ E, assuming that the field A does not depend on 
x~ leads to the conservation of at each vertex coupling the field to the parton. Alternatively, 
this conservation law could be obtained by observing that the transfer of + momentum from the 
plasma to the parton is negligible compared to the large energy E of the parton. 
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potential, where the invariant p + = E component of the momentum plays the role 
of the mass. 

The propagation of the gluon from (t 0 ,x 0 ) to (t,x) is described by the Green’s 
function (x\Q(t,t 0 )\x 0 ) that solves the 2-dinrensional Schrodinger equation (3.19) 
(see e.g. Ref. 108 ). Its solution can be formally represented by a path integral in 
transverse coordinate space: 


(x\g(t,t 0 )\x 0 ) = Q(t-t 0 ) 
where retarded conditions are used and 
U r (t, t 0 ) = Texp 


r 

r e r l , „ , 1 

/ Vr exp 

/ dt r 2 (t ) 

' x 0 

- ^ Jta 


ig f d t'A(t',r(t')) ■ T 

Jtn 


Ur(t,to), (3.20) 


(3.21) 


is a Wilson-line along the path r(t’), with r(to) = and r(i) = x. Between two 
scatterings, the gluon propagation is free, 


(x\Go(t,to)\x Q ) 


0(t - t 0 ) 


27ri(f- to) 6XP 


iE(x — x 0 ) 2 

2(t-t 0 ) 


(3.22) 


The Green’s function g(t,t') plays the role of a unitary time evolution operator for 
a free particle evolving from the position state \x) to the state \x'). In momentum 
representation, 

\p) = J d 2 xe~ zp x \x), / (4j^|P)(Pl = 1 > ( 3 - 23 ) 


we define the amplitude for a gluon of (large) momentum po = {p^ ,Po = 0), present 
in the system at time to, to evolve in the medium into a gluon with momentum 
Vi = (fit > Pi) > where fit = fit = E (to within an irrelevant phase factor), 


M{pi\po) = (pi| g{ti,t 0 ) \p 0 ). 


(3.24) 


The amplitude M. = A4 a b is a matrix in color space, propagating a gluon with color 
b to color a. It is also a (diagonal) matrix in the helicity space. We do not write 
explicitly color and helicity indices to alleviate the notation. 

The free Hamiltonian in 2-dimensional space, reads 

r,2 

(3.25) 




The interaction part of the hamiltonian is the sum of two components, 57~L = bT-L\ + 

SH a , 

(q\6'H 1 \p) = gA(q-p,t), and (q;q'\6H 2 \p) = g(q;q'\V\p). (3.26) 


where 5T~Li represents the elastic scattering of a particle off the potential A(q,t), 
while 57-1.2 represents the splitting of a particle into two. Note that the latter process 
does not conserve the non relativistic energy p 2 /2E, but it conserves the “mass” E. 
When using the propagator (3.20) we effectively treat bHi to all orders. Later, we 
shall calculate the branching probability of a gluon to first order in 5712 ■ Multiple 
branchings will be treated via probabilistic cascades. 
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3.2. Independent multiple scatterings approximation 

The color field of the plasma, through which the jet propagates, fluctuates from time 
to time and place to place. The typical scale characterizing such fluctuations is the 
Debye screening mass tod, which is much smaller than the energy E of the jet. As 
we have seen earlier, the energy E plays the role of the mass in the 2-dimensional 
non relativistic dynamics in the transverse plane. Because this mass is large, the 
motion in the transverse plane is a slow motion, compared to which the dynamics 
of the plasma appears to be very fast. It follows that the interactions mediated by 
exchange of plasma fluctuations are thus nearly instantaneous®. In performing the 
average over the plasma fluctuations, we may then use the following correlator 


(A a (t, x)A* b (t', y)) = nS ab S(t' - i) 7 (* - y) 


(3.27) 


where n is the density of color charges, and 7 (a; — y) will be specified shortly. This 
form of correlator is somewhat analogous to that used in the McLerran-Venugopalan 


model 109,110 that deals with cold nuclear matter of large nuclei. It is also in line 


with the Gyulassy-Wang model, in which the medium is modeled as a collection of 
static scattering centers 28 ’ 29,111 . 

Assuming that it can be reliably calculated from perturbative QCD, the corre¬ 
lator j(x) is given in leading order by 



(3.28) 


which is proportional to the Coulomb elastic cross-section, with q the transverse mo¬ 
mentum transfer. Here we introduced a shorthand notation for transverse momen¬ 
tum integrations, to be used throughout this review: J' = d 2 q/(2n) 2 . A somewhat 
more elaborate version of 7 for a QCD plasma in thermal equilibrium is obtained 
from the substitution 112 q 4 —>• q 2 (q 2 + m 2 D ). 

The picture that emerges form this discussion is one where the jet propagates 
through a random background field, to be averaged over with the correlation func¬ 
tion given by Eq. (3.28). Using another language, the high energy jet collides with 
the plasma constituents at a rate determined by the Coulomb cross-section. Note 
that this depends on the medium solely through the density of color charges, and 
the Debye mass. The motion of the jet in the transverse plane is a slow motion, 
and the interactions with the plasma constituents are effectively instantaneous. 

The treatment of collisions in terms of the propagation in a random field with 
an instantaneous correlation function implicitly assumes that the collisions are in¬ 
dependent. This requirement is satisfied in a weakly coupled quark-gluon plasma: 
The Debye mass is ma ~ gT , with T the temperature. The elastic mean free path 
for soft collisions (with momentum transfer of order mu) can be estimated as 



(3.29) 


g Said differently, because of Lorentz time dilation, the high energy parton has a very poor resolu¬ 
tion in a;+, so that the medium correlations appear to him as effectively local in time. 
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where we have used a ~ < 7 4 /and n ~ T 3 . This is parametrically larger (by 
one power of the coupling constant) than the duration of the interactions, that is 
Ad • 


3.3. Momentum broadening 


The multitude of kicks that it receives from its collisions with the medium con¬ 
stituents affects the propagation of a high energy parton. In particular, it modifies 
its transverse momentum. This leads to the so-called momentum broadening , which 
we now discuss. 

The amplitude (3.24) introduced above allows us to calculate the transverse 
momentum spectrum of a gluon at time t, given the spectrum at time to'. 


d N_ 
d rip 



V(p,t\p 0 ,t 0 ) 


dN 
dn po ’ 


, _ dEd 2 po 

lpo ~ 2E(2 tt) 3 ’ 


(3.30) 


where dN/dfl Po is the spectrum at time t 0 , and V(p,t\po,to) will be referred to 
as the pj_-broadening probability. This is calculated by squaring the amplitude, 
averaging over the initial color and polarization, summing over final polarization 
and color, and finally averaging over the background field. In doing this calculation, 
it is convenient to keep the momentum po in the amplitude distinct from that, po, 
in the complex conjugate amplitude. We get then 

M{p\po) M*(p\po) = N z_ x (Tr(p| Q(ti,t 0 ) |po)(Po| QHh, ti) |p)) 

= {2tt) 2 S(Pq - p) V(p, h |p 0 , to), (3-31) 


where the trace is over color indices and in the last line we have factored out the 
momentum conserving delta function that reflects translational invariance. The 
angular brackets in the first line denote medium average. 

In Eq. (3.31), appears the Fourier transform of a fundamental quantity, which 
we refer to as a 2-ppoint function, (Xi|5'^ 2 ^(t,to)|^o) (see Eq. (3.32)). This is one 
of a collection of n-point functions that play an important role in the discussion. 
This 2-point function, is not a usual field theory 2-point function. It is built from 
the product of a propagator Q of an amplitude and the corresponding propagator 
Q' in the complex conjugate amplitude. This is given, quite generally, by 

(A'i|S' (2) (t,f 0 )|X 0 ) = jv2 1 _ - (Tr (aji|0(ii,t o )|a:o) (x 0 |^ t (t 0 ,ii)|®i)), (3.32) 

where Xi = This two point function describes the evolution from coor¬ 

dinates |X 0 ) = \xq,xq) to \Xi) = \xi,xi) h , of a (fictitious) dipole made of two 
gluons in a color singlet state, one in the amplitude and the other in the complex 

h Some explanation about the notation is needed here. We regard (ti,to) as a matrix in 
transverse space, (Y’|S'( 2 )(£i,£o)|-20j with X = (a?;*), Y = (y\y) with the coordinates x,x,y,y 
attached to Q and Q< as indicated in Eq. (3.32) (the time variables are not considered as matrix 
indices, and S^ 2 \ti,to) is proportional to 6{t\ — to))- 
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conjugate amplitude. From the equation satisfied by the propagator Q, Eq. (3.19), 


one easily deduces the equation satisfied by S( 2 \ or equivalently that satisfied by 
the broadening probability (V{p,t\p n ,t 0 ) = V(p 1 t)) 65 



(3.33) 


where cr(q) is the Fourier transform of the so-called dipole cross-section 

a(x) = 2 g 2 [ 7 ( 0 ) - 7 (a;)], 


(3.34) 


with x corresponding to the transverse size of the dipole. Eq. (3.33) can be solved by 
Fourier transforming to coordinate space. The r.h.s of Eq. (3.33) becomes local in x , 
and the solution exponentiate. Moving back to Fourier space, it is straightforward 
to get (for t = L) 



(3.35) 


which relates the gluon scattering cross-section at high energy with the S-matrix 
for a dipole passing trough the medium, in the independent multiple scattering 
approximation. 113 

In the diffusion limit, namely when the final momentum p is acquired by a large 
number of small momentum transfers, i.e., q <C p, Eq. (3.33) takes the form of a 
Fokker-Planck equation 



(3.36) 


where the (momentum dependent) diffusion coefficient 



(3.37) 


stands for the jet-quenching parameter, 32 estimated here to logarithmic accuracy: 
the lower cut-off of the logarithmic transverse momentum integration is given by the 
Debye screening mass mo- We shall see later that this coefficient receives potentially 
large radiative corrections. Another article in this volume discusses an independent 


set of corrections to the leading-order expression (3.37) (see also Ref. 114 ). 


Finally, assuming q to be roughly constant, one obtains from the diffusion equa¬ 
tion (3.36) that the transverse momentum squared acquired by a high energy gluon 
traversing a length L of the medium is 


(Px)typ = Q L - 


(3.38) 


It can also be shown that, in the harmonic approximation, the dipole cross-section 
(3.34) relates to the quenching parameter as follows 



(3.39) 
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3.4. Quasi-local gluon branching 

We turn now to inelastic processes that are responsible for the formation of the 
gluon cascade. This will be handled in first order in the Hamiltonian 8 H 2 ■ We 
follow here the derivation presented in Ref. 65 

We aim to extract the inelastic transition probability V 2 {k a , kb, L\po, to), that 
relates the 2-gluon spectrum to the single gluon spectrum 


dN 



V2(ka,k b ,L\p 0l t 0 ) 


d N 

dn po 


(3.40) 


This is obtained from the amplitude for a gluon of momentum po = (E,po) at time 
to to propagate until the time t\ where it has acquired the transverse momentum p\ 
and branches into two gluons with momenta q\ = ( zE,qi ) and q[ = ((1 — z)E,q[) 
respectively. The offspring gluons propagate then until they escape the medium at 
L with momenta k a = ( zE , k a ) and kb = ((1 — z)E , kb). This amplitude is given by 
(we omit the color indices to simplify the notations), 


9 

2 E 


dti f (fc a |£(L,fi)|<?i) (fcf,|0(L,ti)|q() {qi; q[\V\pi)(pi\Q(ti,to)\po). 

dqiq'xPi 


(3.41) 


It is represented in the upper diagram of Fig. 3.4. The 3-gluon vertex is given by 


(qa: q'b\V\pc) = 8 (q - q' - p) 2 (q)e^q')e k (p)V jk (q - zp , z), 

(3.42) 


with 

T ijk (Q, z) = —-— Q k 8 ij - Q l S jk + -Q j 5 ik . (3.43) 

and 2 = u a /E is the energy fraction of the parent gluon carried by gluon a. The 
particular combination of momenta Q = q zp, reflects a Gallilean invariance of 
the effective 2-D dynamics. A detailed derivation of this formalism from Feynman 
diagrams in light-cone gauge, A + = 0, is given in Appendix A of Ref. 65 Eq. (3.41) 
provides the order g correction to the evolution operator Q. 

By squaring the amplitude (3.41) in order to obtain the 2-gluon production dis¬ 
tribution, we are left with the average over the gauge field. Because the correlations 
of two A fields is instantaneous, this operation is conveniently visualized by plotting 
the amplitude and its complex conjugate as in Fig. 3.4). By performing the aver¬ 
age over the field configurations with the help of Eq. (3.27), we identify then three 
factors corresponding to the three different time intervals, involving two, three and 
four propagators respectively. The 2-point function, S^ of the first interval [io,ii] 
is related to the broadening probability, cf. Eq (3.30). The 3-point function, S ^ 
of the second interval, will determine the branching rate as will shall see 
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to 


tc\ 


Fig. 1. A diagrammatic representation of the time evolution of the amplitude for a gluon branch¬ 
ing given by Eq. (3.41) (upper diagram) and the complex conjugate amplitude (lower diagram) 
involved in the calculation of the branching probability Eq. (3.45), form to to L. The branching 
takes place at ti in the amplitude and t 2 in the complex conjugate amplitude (see text for further 
explanations). 


shortly. Finally, the last interval, involves a 4-point function, S ^ which factorizes 
into two 2-point function in the limit of a short formation time, 65 

(■ Kk b ; k a k b \S^(L, t 2 )\q 2 q 2 , 92 ^ 2 ) - 

<%2 - q' 2 )S{q 2 - q 2 )V{k a ,L\q 2 ,t 2 )V(k b ,t\q' 2 ,t 2 ), 

(3.44) 

up to terms of the order of t^/L that we neglect. One can indeed show that the 
color correlations between the offspring gluons are exponentially suppressed after a 
time of the order of the branching time 7b r , as they move apart. This property is 
closely related to that of color decoherence that will be discussed in more detail in 
Section 6. 

By integrating over the intermediate states, summing over the final polariza¬ 
tion vectors and averaging over the initial ones using the completeness relation 
Ta=i 2 e \(^) e *x (k) = one finds that the amplitude squared depends on two 
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quantities, the broadening probability and the reduced 3-point function as follows 


V 2 (k a ,k b ,L\p 0 ,t 0 ) = -^2Re [ dh I ' dt 2 [ T i ^(Q 1 )P^(Q 2 ) 

Jo Jo JQi ,Q 2 .1 

V(k a , t\q 2 , t 2 )V{k bl t\q' 2 ,t 2 ) S (3) ( Q 2 , QiZ; ii, t 2 ) V{pi,t\po, t 2 ), 

(3.45) 


where the 3-point function £o 3 ) reads (cf. Ref. 78 ) 

S (3) (Qi,Q 2 ,l;h,t 2 ) = J du 1 du 1 dve lUlQl - lU2Q2 ~ l ' vl 


Vu 


x exp 




dt [cr(u) + cr{v — zu) + cr(v + (1 


-)«)] j . 

(3.46) 


with Q i = qi — zpi, Q 2 = q 2 — zp 2 and l = p 2 — p±. Momentum conservation at 
each of the two vertices implies that pi = q± + q[ and p 2 = q^ + q 2 - The function 
S( 3) is an important quantity that encodes the evolution, from t\ to t 2 , of the the 
system of two gluons with energies zE and (1 — z)E respectively, and one gluon with 
energy E. The transverse vector u corresponds to the transverse distance between 
the offspring gluons a and b, while v corresponds to the coordinate of the center of 
mass of the system. 

Furthermore, the sum of gluon momenta in the amplitude and the complex 
conjugate amplitude vanishes at any given time, <72 + q 2 = q 2 + q 2 , Pi = P2 and at 
each vertex we have q[ = Pi — qi and q' 2 = p 2 — q 2 . The contraction of the three 
gluon vertices in the amplitude and complex conjugate amplitude yields 

ir^(Q 1 )r i ^(Q 2 ) = -tt~ \P( z )(Qi ■ Q 2 ), (3.47) 

2 Zy 1 Z) 


where 


P{z) = N c 




(3.48) 


is the gluon-gluon Altarelli-Parisi splitting function. 95 Eq. (3.45) can be further 
simplified in the short formation time approximation, i.e, Tb r <C L. 

The difference between t\ in the amplitude and t 2 in the complex conjugate 
amplitude, r = At = t 2 — t\ is at most of the order of Tb r - The 3-point function 
decays exponentially for r > Tb r - It follows that the upper bound of the integral 
over r can be send to infinity, 

pL pt2 rL pt2 pL nOO 

/ dt 2 / dti= dt 2 / dr« 1 dt 2 dr. (3.49) 

Jo Jo Jo Jo Jo Jo 

Similarly, the relative transverse momenta generated in the branching, Q 1 , Q 2 , l 
are typically of the order of fcbr ~ VQ T br and thus, can be neglected compared to 
the external momenta. This is the collinear branchings approximation. Hence, one 
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can integrate over these internal transverse momenta and r, and this allows us to 
define the transverse momentum integrated kernel, 33,115 1 


JC(z,E) = 


P(z) 


[z(l-z)E]* 


Re 


/dr/ (Qi-Q2)S^(Q 1 ,Q 2 ,l,t 1 ,r), 
Jo JQ 1 Q 2 I 


(3.50) 


In the harmonic approximation, we can ignore the logarithmic dependence of q that 
appears in the dipole cross-section (3.39), and the reduced 3-point function (3.46) 
can be explicitly calculated. This yields for the kernel (3.50), 


JC(z,E) 


[1 -z + z 2 ] 5 / 2 U 

[^(1-z)] 3 / 2 \l E' 


(3.51) 


Finally, the transition probability (3.45) takes the compact form 65 (letting t\ = t ) 


V 2 {ka,kb,z]L\p 0 ,to) 

= 2g 2 z(l-z) / dt'P(k a ,L\zp 1 t)V(k b ,L\(l 
Jo 


z)p,t) JC{z,E)V(p,t\po,t 0 ). 

(3.52) 


When the branching time is comparable to the size of the system, Tb r ~ L , finite 
size effects are no longer negligible but in this regime only a single branching oc¬ 
curs. 65,69 The corresponding branching probability integrated over transverse mo¬ 
mentum generalizes Eq. (2.6) including finite energy corrections, that is z ~ 1 — z, 
can be found in Refs. 33,115,116 Also, the generalization to arbitrary parton species 
is straightforward and can be recovered from the aforementioned references. 

This is a leading order result that describes the instantaneous branching of a 
hard gluon inside a medium of size L. The probability distributions (3.35) and 
(3.52) can serve as building blocks for the in-medium gluon cascade that we shall 
present in Section 3.5. 


3.5. Independent multiple branchings approximation 

The probability distribution (3.52) is parametrically of the form 

V 2 ~ —L, (3.53) 

Tbr 

where d/rb r is the rate of (quasi-instantaneous) branching, and the factor L reflects 
the fact that the branching can occur anywhere along the medium, i.e., 0 < t < 
L. Consider now the branching probability Vs, where two successive branchings 
produce three gluons in the final state. The two branchings occur successively at 
times t\ and t 2 in the amplitude and at times t[ and t ' 2 in the complex conjugate 
amplitude. When the time intervals [ti,^] and [£ 2 ,^ 2 ] d° n °t overlap, that is, 
t\ < t,\ < t 2 < f(, • the resulting successive integrations over the different times 

~ /q\ 

1 The vacuum part of the 3-point point function, S q , is implicitly subtracted. 
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yields Lr r^ r , provided Tb r ~ ribr ~ T 2 br "C L. In this case, the branchings are 
independent from one another and the 3-gluon probability factorizes as 

P” dep ~ (P 2 ) 2 ~ (a^j ■ (3.54) 

On the other hand, interferences between the subsequent branchings may occur 
when the time intervals [ti, t \] and [^ 2 ,^ 2 ] overlap. For example this may happen 
when t± < t 2 < t'[ < t ' 2 with ti,T 2 < Tb r L. The corresponding time integrals 
yield a factor L r^ r . In that case, the contribution to P 3 from the interferences, 
P" lter , is suppressed by a factor t\, t /L with respect to the factorized form (3.54), 
and it is obvious to generalize the argument to all possible time configurations. 

For consistency, and to complete the picture of independent multiple branching 
approximation, it is of course necessary that the time between successive branch¬ 
ing be large compare to the branching time. But this is guaranteed for inclusive 
quantities at weak coupling since t*(w) = Tb r (w)/d Tb r (w). 

It follows that, as long as the length L of the system is much larger that the 
branching time Tb r , successive branchings of a gluons in a large medium can be 
treated in first approximation as a probabilistic cascade. The properties of this 
cascade can be deduced from a set of n-point functions for which a generating 
functional was constructed in .' 8 

The two building blocks needed in this construction are the elastic and collinear 
branching rates, or equivalently the probabilities Pi and P 2 given by 

Pi{k\p) = 2(2 n)E 6 (u - E)V(k , t\p, t ), 

P 2 (ki,k 2 \p) = 2(2ir)E5(u > 1 + to 2 - E)S(p-k 1 - k 2 )5(ki - zp)V 2 [ki, k 2 , z; t\p, t 0 ). 

(3.55) 



Fig. 2. Graphical representation of the elastic (left panel) and inelastic (right panel) rates. For 
collinear branching, q = k/z and 1 = 0. 


In this review, we shall focus solely on the lowest of the aforementioned n-point 
functions, namely the single-inclusive gluon distribution 


D(u, k) = oj 


d N 


dojd 2 k 2(27t) 3 


1 


Y. / II d flinP n (k,ku--- ,k n -i). (3.56) 


i= 1 
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It obeys the evolution equation 

d 
dt 


N n f 

D(x,k,t) = —J— / a(l,t)D (x,k — l,t) 

2 Jq 

/ x \ ^(x k \ 


a s d z 
Jo 




— JC (z, xE) D (x, fc, t) 


(3.57) 

The first and second terms in the r.h.s. correspond, respectively, to the elastic and 
the branching rates represented in Fig. 2, with in the latter case 1=0. 


3.6. The gluon distribution beyond the collinear approximation 

In constructing the probabilistic picture in the previous section, we have explicitly 
assumed that the dominant contribution to the transverse momentum broadening 
along the cascade comes from elastic scatterings. Some increase in transverse mo¬ 
mentum can also arise from the branching processes. 

Let us then estimate the contribution of multiple radiations to transverse mo¬ 
mentum broadening. Recall that the typical transverse momentum generated by 
one radiation of frequency w is fcb r (w) ~ (wg) 1 / 4 . Combining this with the radiative 
spectrum, ujdl/duj ~ L/t*(pS) (cf. Eq. (2.5)), we obtain the following parametric 
estimate for the total squared transverse momentum generated in splittings 

( fc I)typ,rad~ I ^T k br M ~ aqL In J-. (3.58) 

Ju BH W A el 

where we have integrated over the relevant range of gluon energies from ubh to 
u) c , and we have used the relation wbh = <?A e i. Compared to the typical transverse 
momentum acquired by elastic scattering, the latter estimate is suppressed by a 
factor a, yet enhanced by a potentially large logarithm In L/A e i. We will see that 
this not the dominant contribution. A more detailed analysis reveals a double 
logarithmic enhancement' 5 that is missed by the estimate we just did. We shall 
return in the next section to a detailed discussion of the origin of this correction. 
Here, we just indicate how it can be simply obtained by taking into account the 
momentum broadening taking place during the branching process. This forces us 
to go beyond the collinear approximation, and hence perform a calculation whose 
validity is not immediately obvious, but will be justified later. 

We note that the 3-point function defined earlier contains the information about 
the transverse momentum broadening accompanying a single branching. Hence, in¬ 
stead of integrating over l and Q, as we did for Eq. (3.50), we define the unintegrated 
kernel 

K(Q,l,z,E) = P{Z \ Re r dr [ {Q 1 ■ Q)S^(Q 1 ,Q,l,t 1> r). 

[z( 1 - z)E ] 2 J o J Q 


(3.59) 
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One can thus generalize Eq. (3.60), 

^D(x,k,t) = a s J d zj ^ K (q,1,z,^eJ D 


— K. ( Q , Z, z, xE ) D (x,k — Z, t) 


N r n 


j a(l)D (x,k — l,t) . (3.60) 


where Q = k — z{q +Z). We shall now follow the strategy that we adopted in order 
to reduce Eq. (3.33) to the diffusion equation (3.36). This involves an expansion 
around the large momentum k of the followed gluon. The momenta Q and Z, which 
are at most of the order of fcb r = z( 1 — z)qE , are small compared to k which is 
typically of the order of qL. The expansion of the gluon distributions around k 
yields for the first term of Eq. (3.60) 


D 


x k — 6k\ 

) 


= D 



-6k 




1 

2 ! 


d d 

5k ' 5k3 ww D 

r\j'i r\j j 



H- 

(3.61) 


where we have set 6 k = Q + zl. One expands similarly D {x,k — V). It is easy 
to see that the leading terms reproduce Eq. (3.57). The linear terms vanish upon 
angular integration. Remain the quadratic terms, whose contribution can be cast 
in the form of the diffusion term, thereby exhibiting a correction A q to the jet 
quenching parameter. For consistency, we shall also simplify the collision term by 
using the diffusion approximation. At this point, we anticipate that evaluation of 
the correction A q meets with logarithmic divergences. These arise from the region 
2 ~ 1. To the leading-logarithmic accuracy, we can set 2 = 1 everywhere, except in 
the dominant singularity. The dominant contribution to A q can be then written as 

A q{k 2 ) = 2a s [ dz [ [(Q + l) 2 -l 2 ])C{Q,l,z,xE), (3.62) 

Jx JQ,l 

where the k 2 dependence arises from the integration boundary Q 2 ,l 2 kr. The 
complete calculation of the integral (3.62) is presented in the next section. One 
gets 78 


A q(k 2 ) 


a s N c 

2 n 



(3.63) 


where 1 is the maximum energy that can be extracted from the medium in a single 
scattering (e.g. r 0 = 1/T for a weakly coupled plasma with temperature T). We 
return to this result in the next section. Finally, Eq. (3.60) can be recast in a form 
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that is similar to that of Eq. (3.57) (after performing the diffusion approximation), 78 

^D(x, k, t) = ^ (Jl^j [(g(fc 2 ) + A q{k' 2 ))D (. x, k - q, t)] 

+ a s J dz — K, (z, -Pj D —, t\ - 1C (z, xE) D (x, k, t) 


(3.64) 


where p ±-broadening due to soft gluon radiation is now taken into account effec¬ 
tively via a redefinition of the quenching parameter q —> q + A q. 

We are now ready to revisit the standard estimate of the typical momentum 
broadening Eq. (3.38), making use of the redefinition of the quenching parameter 
above, substituting k 2 ~ qL in Eq. (3.63), we obtain 


(fcl)typ - fa 



ot s N c 2 L 
-In — 

7T Tq 


(3.65) 


This result agrees with that obtained in Ref. 75 using a different approach^ Note 
that we obtain a double logarithm of the medium length, as opposed to the heuristic 
estimate (3.58), which involve a single logarithm. Although, the soft (short time) 
logarithmic enhancement is correctly captured in this estimate, the single scatter¬ 
ing that contributes to the branching probability that in turn leads to the second 
logarithm, is missing in the BDMPS-Z spectrum that is used Eq. (3.58). 


4. Renormalization of the jet quenching parameter 

The correction to the jet quenching parameter discussed in the previous section 
points to the existence of potentially important radiative corrections. The calcu¬ 
lation that we just presented focussed on a correction that is singular, and this is 
what allowed us to retain in the branching kernel contributions that should not be 
kept a priori since, if they were not divergent, they would be of the same order of 
magnitude as terms that have been consistently neglected in arriving at Eq. (3.57). 
In this section we discuss a more systematic approach to the calculation of radiative 
corrections. Still, at present, the program of calculating the radiative corrections is 
incomplete. 74,75 ’ 78,79 ’ 117 ’ 118 As we shall see, one is able so far to complete the calcu¬ 
lation only for the leading singular part. Furthermore, it is only this singular part 
that can be interpreted as a correction to the jet quenching parameter; less singular 
corrections would presumably be non local, thereby spoiling their interpretation in 
terms of a correction to the jet quenching parameter. 

J Note that NLO corrections to p ±-broadening have been also investigated in Deep Inelastic Scat¬ 
tering in the High-Twist approach. 77 
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4.1. Radiative corrections to momentum broadening 

We start by revisiting the calculation of the momentum broadening. Recall that this 
is described by the 2-point function (ti,to), Eq. (3.32), for which the following 
equation is easily established 

S {2 \t i, t 0 ) = So 2 \h,t 0 ) + j dt 3 J dt 2 S^\t ll t 3 ) E (2) (f 3 ,t 2 ) S^\t 2 ,t 0 ). (4.66) 

In this equation, the instantaneous interaction kernel E^ 2 ) is related to the dipole 
cross-section (3.34), 

(a;|E {2) (t 3 ,t 2 )|y) = - y)S(t 3 - t 2 )a{x). (4.67) 



to to ^3 t\ in to £3 1 i 

Fig. 3. Real (left) and virtual (right) contributions to the radiative correction to the 2-point 
function. There are other diagrams that are not shown, corresponding to a different time ordering 
of the emission and absorption vertices in the real term, and a diagram where the virtual correction 
is attached to the lower line. The transverse vector v stands for the transverse size of the frozen 
hard dipole (represented by the straight lines), whose constituents carry an energy E , and u 
corresponds to the time-dependent transverse distance between the radiated gluon (represented 
by a wavy line) and its radiator whose dynamics is encoded in the 3-point function The 

integration over u ~ q~ 1 in the region u^> v yields the second logarithmic enhancement of the 
radiative correction to q (cf. Eq. (4.74)). 


We are interested in the radiative correction to E*- 2 ). The corresponding dia¬ 
grams are displayed in Fig. 3. The important par is that involving the three gluons 
evolving between the times t 2 and t 3 , and to which we refer as a 3-point function. 
The dominant, double logarithmic contribution corresponds to the case where the 
radiated gluon (the wavy lines in Fig. 3) is soft, i.e., its energy w is « C F. A 
simplification then occurs in the calculation of the diagrams since we may assume 
that the transverse coordinates of the hard gluons (in the amplitude and its com¬ 
plex conjugate) are frozen during the lifetime of the fluctuation, i.e., between t 2 (in 
the amplitude) and t 3 (in the complex conjugate amplitude). The calculation leads 
then, for the real correction, and using a symbolic matrix notation, 

AS^(h,to) = - r dt 3 [ 3 dhSWih^AXW^t^SWfato), (4.68) 

Jto Jto 
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with 


(X 3 \AY 2 \t 3 ,t 2 )\X 2 ) 

= 6{X 3 - X 2 ) a PY 2 Re J ^ d U2 ■ d U3 S {3) (u 2 ,u 3 ,v;T)\ U2=0iU3 ^ v , 

(4.69) 

where r = t 3 — t 2 . The 3-point function S < ' 3 \u 2 ,u 3 ,v,t) encodes the dynamics 
in the transverse plane of the system made of a frozen dipole and a soft gluon. It 
corresponds to the soft limit, i.e. z = to/E <C 1, of Eq. (3.46). Here, v stands for 
frozen transverse size of the hard dipole during the time interval r, while u 2 = 0 
and u 3 = —v stand for the transverse separation between the radiated gluon and its 
radiator (the hard gluon in the amplitude) at the times t 2 and t 3 respectively (see 
left panel of Fig. (3) for an illustration). The virtual correction is given by the same 
formula in which u 3 = 0 since the gluon is emitted and reabsorbed by the same line. 
Anticipating on an argument that will be provided later, we treat the correction 
^ 2 \t 3 ,h) effectively as a local correction, i.e., set formally T,( 2 \t. 3 ,t 2 ) a S(t 3 —t 2 ). 
After replacing in Eq. (4.68) the factors external to the radiative correction 
by free 2-point functions (which are independent of time), and performing the t 2 
integration, we get 

(X 1 \ASW(t 1 ,to)\X 0 ) = -(t 1 -t 0 )6(X 1 -X 0 )^L J^J d TK(v,r), 

(4.70) 


with 


K(v,t) = 2Re f V - e U v ) (q 2 ■ q 3 )S^(q 2 ,q 3 } l-r), (4.71) 

Jq 2 ,q34 v 7 

where we have performed a Fourier transform of the 3-point function. Besides an 
integration over r, we recognize in Eq. (4.71) the soft limit, i.e. z <C 1 of the 
kernel Eq. (3.59). Note the vanishing of this expression when v —> 0. It results 
from a cancellation between real and virtual terms, reflecting the property of color 
transparency of a dipole cross section. 

The identification of the correction to the dipole cross section proceeds by com¬ 
parison with the integral equation (4.66) (more properly, its generalization, whose 
validity relies on the locality assumption) for fixing the relation between Act and 
AE( 2 ). One then gets, leaving the bounds on the r integration unspecified for the 
moment, 


N c n 


A a(v) 



dw 


dr K(y, t) 


1 

4 


v 2 Aq, 


2 


(4.72) 
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where the last equality stems from Eq. (3.39) expressing A a(v) at small v in the 
harmonic approximation, by expanding the phases in Eq. (4.71). We use this re¬ 
lation to interpret the correction Acr(v) as a correction to the parameter q. Note 
that in this approximation the resulting integrations over the transverse momenta 
q 2 and l are bounded from above by v~ 2 . 

The explicit calculation of A q, using the expression (3.46) of the reduced 3-point 
function, yields' 8,79 


A q 


:N r 


2Re 



dT 


in 3 

sinh(Qr) 


a s N c f dw f n dr „ 
— ~ Q ’ 



4 

sinh 2 (f It) 


(4.73) 


which as anticipated exhibits a double logarithmic divergence, when r —► 0 and to —* 
0. As clear from Eq. (4.73) the r integral is bounded at the upper end by ~ 
Tbr(w) = i Jw/q (cf. Eq. (2.7)) corresponding to the onset of the multiple scattering 
regime and the LPM effect: the relevant gluon fluctuation experiences a single 
scattering with the medium constituents. In order to systematically account for the 
boundaries of the double integral, k it is in fact more convenient to change variables, 
from (w, r) to (q, r), with r = w/q 2 the formation time of the radiated gluon, and q 
its transverse momentum which can run up to p 2 = v~ 2 (the logarithmic transverse 
momentum integration can be rephrased in terms of the transverse coordinate of 
the radiated gluon q -1 ~ u v, see Fig. 3). We obtain then 



where we have explicitly indicated the scale dependence of q. The boundary cor¬ 
responding to the region of multiple scattering now appears as the lower bound of 
the q integration, q 2 fc 2 )r = qr. Since medium-induced gluons forms inside the 
medium, the largest value for r is the length of the medium L ~ t — to. As for the 
lowest value to, it can be interpreted as the inverse of the largest energy that can 
be extracted from the medium through a single scattering. 

For a constant q in the integral, one can easily perform the integrations, and by 
keeping the leading contributions, we recover the result first derived in,' 5, ' s 

(4J5) 

where we have used the fact that p 2 ~ qL. 


k For a detailed discussion on the boundaries of the logarithmic integrals see Ref.' 5,118 
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4 . 2 . Locality of radiative corrections 

So far we have discussed a single radiative correction. The logarithmic phase space 
for radiative corrections extends from r 0 to L which makes them explicitly non-local. 
One may therefore question why we have argued that these corrections could be ab¬ 
sorbed into a redefinition of q , considered to be a local transport coefficient. Indeed, 
the quenching parameter makes sense strictly speaking only in the approximation 
of independent multiple scatterings. This requires the duration of each interactions 
to be much smaller than the mean-free path. 

The key lies in the nature of the logarithmic integral whose boundaries do not 
affect the overall multiplicative constant. Thus, to the extent that one restricts 
oneself to a leading order calculation, one can proceed as if the lifetimes of the 
fluctuations involved in the radiative corrections were small, and treat them as 
local. The coefficient of the leading double logarithm is calculated correctly, and 
corrections coming from overlapping contributions will be subleading. In short, 
the leading corrections can be treated as being effectively local in time and thus 
independent. ' 9 This argument is what allows us to generalize Eq. (4.66) to a E(f, t') 
that includes the radiative corrections, to treat the double logarithmic correction to 
£(f,f') as effectively local, and multiple radiative corrections as independent. As we 
have seen, such corrections can be interpreted, in the coordinate space description, 
as a modification of the dipole cross-section o(v) —► a(v) + Aa(v, r max ), where 
T max = t\ — to, or in momentum space as a modification of the equation for the 
momentum broadening probability' 9 



(4.76) 


which generalizes Eq. (3.36). In both cases, the radiative corrections are accounted 
for by a correction to the jet quenching parameter. 

4 . 3 . Radiative corrections to Energy Loss and Universality 

We now address the question of whether the absorption of the dominant logarithmic 
contributions to the radiative corrections in a redefinition of q is universal, that is, 
whether it holds for other observables than the broadening probability. This has 
important implications as it would affect the branching rate in Eq. (3.64), that 
controls the number of final state gluons. In order to address this issue, we consider 
radiative corrections to the BDMPS-Z spectrum (2.5) that can be recovered in limit 
w = zE <C E of the branching rate (3.50). 

The BDMPS-Z spectrum (2.5) of radiated gluons with frequency w -C E, can 
be calculated from the reduced 3-point function (3.46), 33 according to 1 



(4.77) 


'The vacuum part of the 3-point point function, sjf 1 , is implicitly subtracted. 
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Here the variable t runs up to ~ L. This spectrum is valid in the large medium 
length limit where the gluon branching time r^ r = ^ui/q <C L and the integration 
over the gluon formation time t < L is suppressed exponentially beyond Tb r . The 
transverse coordinates u and u' correspond to initial and final coordinates of the 
radiated gluon, and v corresponds to the size of the energetic dipole that radiates 
the gluon (see Fig. 4). This is why we can integrate r up to infinity in Eq. (4.77). 



tn tn + T 

Fig. 4. Diagrammatic illustration of the BDMPS-Z spectrum (4.77) that corresponds to the 
radiation of a soft gluon with energy u. We denote by u the transverse position of the radiated 
gluon with respect its radiator, the energetic dipole whose size is denoted by v. 


The 3-point function S^ is given explicitly by Eq. (3.46) in the harmonic ap¬ 
proximation, and t — t\ — to, with t 0 and t\ denoting the times of the emission in 
the amplitude and its complex conjugate, respectively (see Fig. 4). The BDMPS-Z 
spectrum of Eq. (2.5) is easily recovered from this expression by performing the in¬ 
tegrations over the transverse momenta and over r. 65,73 As is clear from Eq. (3.46) 
the reduced 3-point function depends explicitly on q. One can show that the lead¬ 
ing radiative corrections do not modify S except for a change in the value of the 
parameter q, the correction to q being, besides, the same as that calculated in the 
previous section for momentum broadening. 

Quite generally, we are concerned with the radiative corrections of the 3-point 
function whose graphical interpretation is given in Fig. 4. The diagram dis¬ 
played there corresponds typically to a branching process where a gluon with initial 
energy E (represented by the lower two thick lines in the amplitude and the conju¬ 
gate amplitude, respectively), emits a gluon at to in the amplitude and reabsorbs 
it in the complex conjugate amplitude at t\ = to + t, with energy oj = zE. The 
correction to the 3-point function (in the limit v — > 0 as required by the integration 
over l in Eq. (3.46)) can now be written by analogy with Eq. (4.68) as (see Fig. 5 
for an illustration) 


A5' (3) (t 1 ,t 0 ) - - / 1 dt 3 [ 3 A^ 3 \t 3 ,t 2 ) S^ 3 \t 2 ,to), (4.78) 

Jt 0 J t 0 
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to tv tv + T 7 to + t 

Fig. 5. Diagrammatic illustration of the radiative correction to the BDMPS-Z spectrum (4.77). 
The gluon c Y is integrated over and account for real and virtual contributions. Here one topology 
among 9 is displayed corresponding to the radiation of the gluon u' from gluon lj and its absorption 
by the gluon E in the complex conjugate amplitude. The labeling is as in the previous figure, with 
here, in addition, r denoting the transverse coordinate of the gluon fluctuation a/. The integration 
over r ~ q~ 1 in the region r u yields the second logarithmic enhancement of the radiative 
correction to q (cf. Eq. (4.73)). 


where 


(X 3 \A^ 3 \t 3 ,t 2 )\X 2 ) = 6(X 3 -X 2 ) a ^ J ^ K(u,t '). 

(4.79) 


At this point we proceed as for the 2-point function, and treat the correction 
AI]( 3 )(f 3 , £ 2 ) as a local correction. By comparing with Eq. (3.46) in the limits 
z —> 0 and v —>• 0, one then gets 


N c n 

2 


Aa(u) 



d J 

773 


K(u, t') , 


(4.80) 


Thus, to double logarithmic accuracy, the radiative corrections to the reduced 
3-point function are accounted for by correcting the dipole cross-section and thus 
the jet quenching parameter. This result applies to the particular 3-point function 
involved in the calculation of the gluon spectrum (4.77). In order to evaluate the 
corrected spectrum, we need to perform the integration over r in Eq. (4.77), replac¬ 
ing q by q + Ag(r max )- However, since the correction to q is computed to double 
logarithmic accuracy one can simply replace the variable r max in Eq. (4.80) (see also 
Eq. (4.74)) by its typical value in the radiation process, i.e., r max ~ Tb r = \Jui/q. 
The integral over r in the BDMPS-Z spectrum (4.77) can then be performed as for 
the case with no radiative correction. Doing so, we obtain lfl 


dJ 

dwdf 


S JV, 


A q 


u> 


7T 


(4.81) 
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where for a constant g one gets, from Eq. (4.74) letting p 2 ~ fc^ r (w) = y/cJq, 



(4.82) 


It can be shown that this result extends to the full kernel (3.50) in the large N c 
limit. 79 These results were recently confirmed by an alternative calculation of the 
radiative corrections to the mean-energy loss to double logarithmic accuracy. 11 ' 

4 . 4 . Renormalization of the quenching parameter 


For large media, as soon as aln 2 (L/r 0 ) ~ 1, one has to resum the double logarith¬ 


mic power corrections. Unlike the previous resummation of independent multiple 
radiative corrections, this now involves radiative corrections that are correlated to 
each other. To understand how this resummation proceeds, we denote the standard 
leading order definition of the jet-quenching parameter by q 0 and we note that the 
first correction to the jet-quenching parameter, qi(r, k 2 ) = Ag(r, k 2 ) is proportional 
to the 3-point function, which is itself a function of the leading order g 0 . As 

we have shown, under radiative corrections, the 3-point function gets renormalized 
by a simple modification of g, that is, S’® [go] —t S^^go + gi]- This allows us to 
compute the second correction from Eq. (4.74), 



(4.83) 


One sees emerging a self-similarity that results from the separation of time scales 
involved in the computation of the leading logarithms. The structure of the first 
double logarithmic corrections being set, the next corrections that yield double 
logarithms will follow the same systematics, with successive gluonic fluctuations 
ordered in formation time ^ <C q < ... < r„ s T max , or in transverse size ro S> 
ri ... r n = r max , or in transverse momentum mn qi •C ••• q n = k. 
A diagrammatic illustration is given in Fig. 4.4. The difference with the standard 
Double-Logarithmic Approximation (DLA) is the limits of logarithmic phase-space 
set by the LPM effect since, i.e., multiple-scatterings since in the DLA only a single 
scattering contributes, which imposes that the formation time of a fluctuation to be 
smaller than the BDMPS-Z formation time, or in terms of our transverse momentum 
variables, q 2 q$T m . The following equation resums the double logarithmic 

corrections to all orders 



(4.84) 


with some initial condition g(ro, k). We have let the coupling running at the trans¬ 
verse scale q. The important feature of this equation is that it predicts the evolution 


m The resummation of double logarithms in the quenching parameter was postulated earlier in 
Ref. 119 where the LPM suppression was not taken into account. 
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of the jet-quenching parameter from an initial condition q 0 (which can be computed 
e.g. on the lattice, or to leading order in a 3 , which implies, <z(t 0 ) = q 0 as given 
by the leading order result (3.37). The r 0 cut-off that was introduced to cut the 
logarithmic divergence in the radiative corrections, can be seen as a factorization 
scale. The solution to this equation was derived in 75 for the pt-broadening in the 



Fig. 6. A diagrammatic representation of the resummation of double logarithmic contributions 
for a single scattering. Here three gluon emissions are depicted. In the double logarithmic approxi¬ 
mation successive radiations are strongly ordered in formation time, i.e. to <C t\ <C T2 <C T3 L , 
and in transverse sizes, i.e., 1/mo ~ to r\ 7*2 T3 v. For the original dipole the 
subsequent radiations look as if they were instantaneous. 


case where qo = 9 ( 70 ) is constant and for a final r = L and k 2 = q$L , merging the 
2 independent variables at the end of the evolution. The solution reads 

q(L) = -b/i ^2\/dln q(r 0 ). (4.85) 

V a V r o / 

A semi-analytic analysis of this equation including the running of the coupling 
has been performed in Ref. 120 For large L, the quenching parameter scales like 
q(L) ~ L 7 , with the anomalous 

7 = 2-\/d . (4.86) 

Interestingly, the resummation of large double logarithms modifies the scaling of 
the energy loss with L , 

A E ~ i 2+7 , (4.87) 

a correction that seems to fall between the standard small coupling result, A E ~ L 2 
and the strong coupling result obtained with the help of the AdS-CFT correspon¬ 
dence in J\f = 4 SYM theory , 121,122 A E ~ L 3 . 

Finally, let us make now a rough estimate of the renormalized quenching param¬ 
eter. The standard perturbative estimate yields a value of about go ~ 1 GeV 2 /fm . 32 
For a s ~ 0.5, L ~ 5 fm and T = Tg " 1 ~ 500 MeV, Eq. (4.85) yields a sizable increase 
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by a factor 2, q ~ 2 GeV 2 /fm. This result is in the ballpark of the jet quenching 
parameter values extracted from the data. 123 


5. Energy flow and angular structure of the in-medium cascade 


We return now to the medium induced cascade, and investigate some of its char¬ 
acteristic features as revealed by the analysis of the inclusive distribution function 
D(x, k ). We shall discuss in particular the energy and angular distributions of the 
emitted gluons. As we shall see, these properties make the BDMPS-Z cascade very 
different from the more common parton cascades in vacuum that are described by 
the DGLAP equation. 

Our starting point is Eq. (3.57) for the inclusive distribution. Since the trans¬ 
verse momentum k remains small compared to the longitudinal momentum, it is 
actually convenient to use angle variables rather than transverse momenta. Accord¬ 
ingly, we set 


D(x, 0) = (27t) 2 x 


dlV 
dxd 2 6 


where 


w 

X= E ’ 


and 


uj 


k 

xE 


(5.88) 


(5.89) 


Note that 6 is a 2-dinrensional vector collinear to k , whose (small) magnitude 
equals the polar angle of the emitted gluon with respect to the initial direction of 
the leading particle. In this new variable, Eq. (3.60) reads 


d 

dt 


D(x,0,t) 


with (see Eq. (3.34)) 


d zK{z) 


D(x/z,0,t) D(x,8,t ) 

t*{x/z) t*(x) 



d 2 0' 

(2?r) 2 


a ( O ', x)D(x, 9 — O' ,t ), 


cr(x,0) = ( xE) 2 a(q ) = 


327r 2 a 2 

(xE) 2 


1 



(5.90) 


(5.91) 


An illustration of the gluon cascade described by Eq. (5.90) is given Fig. 5. Note 
that the locality (in angle) of the splitting term reflects the effective collinearity of 
the splitting. 

In addition to changing to angular variables, in Eq. (5.90) we have explicitly 
factorized the kernel (3.50) into 


where now 


2 a s JC(z, xE) 


£(*) 

t*(x)’ 


IC(z) 


z 3 / 2 (l-z) 3 / 2 ’ 


and t* (x) 



(5.92) 


(5.93) 
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We shall refer to IC(z) as the reduced kernel. The quantity 


t*{E) = 



(5.94) 


is a characteristic time scale of the BDMPS-Z cascade. After a time of order t*(E), 
most of the initial energy has been radiated into soft gluons (provided L > t*(E), of 
course), hence the often used denomination of stopping time , or stopping distance. 124 
In the following, we often denote t-*(E) simply by f*, except when confusion may 
arise. Note that t*(x) = t*y/x. 

As such the kernel (5.93) does not account for finite size effects that become 
important when i*(x) > L , corresponding to x c > 1 (w c > E) (cf. Eq. (2.6)). In 
the following discussion we neglect these corrections, we thus implicitly restrict jet 
energies to be E ~ w c or smaller (see Refs.' 0,71 for discussions on this matter). 

In the next subsection, we analyze the energy distribution, obtained by integrat¬ 
ing D{x, 0, t) over the angle. That is, we define 

D{x,t) = J -^^D(x,6,t), (5.95) 


which obeys the following equation 44,45 > 66 > 12 5 

9 n( , [ a rt \ \ D ( X / Z > t ) 

— D(x,t) = / d zK.(z) 


D (x, t ) 


t*(x/z) 


t* (x) 


= AD). 


(5.96) 


The angular distribution of the radiated gluons will be discussed next. 


5.1. Energy distribution and wave turbulence 

Since initially the total energy is carried by a single gluon, the initial condition 
for Eq. (5.96) is simply Dq(x) = 5(1 — x). In the regime where the length of the 
medium is so short that at most a single branching can be induced, L < i*(x(l —x)), 
Eq. (5.96) can be solved by iterations, with the first one given by 

Di(x,L) = —x/C(x), (5.97) 

which coincides with the BDMPS-Z spectrum (2.5) at small x. In the opposite 
regime where L 3> t*(x), which corresponds to gluon energies ui <C lo s = a 2 qL 2 , 
multiple branchings are important. A non perturbative solution becomes manda¬ 
tory. 

Eq. (5.96) can be solved exactly for a simplified version of the reduced kernel 
(5.93) in which one neglects the z dependence of the numerator, that is 

zVHi-zy/ 2 ' (5 ‘ 98) 

This turns out to be an excellent approximation. In fact, we shall argue later that 
the exact form of the kernel plays a minor role in the determination of the general 
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t 0 Ati A t 2 A ts A U t 

Fig. 7. Illustration of a gluon cascade that is initiated by a gluon with energy E. Four generations 
are displayed. The branching time Ati ~ t*(a?i), that corresponds to the lifetime of generation 
i, decreases after each branching as in a BDMPS-Z cascade. The inclusive distribution D(x,0) 
measures the probability to find in the cascade, at time i, a gluon with energy xE. The rate 
equation (5.90) describes how this distribution evolves with time t. 


features of the cascade. The solution of Eq. (5.96) for the simplified kernel reads 66,72 

- 2 N L 


D(x,t) = 


y/x (1 — X) 3 / 2 


exp — 7 r 


1 — x 


with r = —. 

1 3k 


(5.99) 


This solution exhibits two remarkable features: a peak near x = 1 associated with 
the leading particle, and a scaling behavior in \/y/x at small x where the x depen¬ 
dence factorizes from the time dependence, i.e. 


D{x, t ) 


(5.100) 


An illustration of this solution is given in Fig. 11, left panel. The energy of the 
leading particle, initially concentrated in the peak at x < 1, gradually disappears 
into radiated soft gluons, and after a time t ~ i* (i.e. r ~ l/\/7r ~ 0.5) most of 
the energy is to be found in the form of radiated soft (a; < 0.1) gluons. This is 
also the time at which the peak corresponding to the leading particle disappears 
(see Fig. 11). At the same time the occupation of the small x modes increases 
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(linearly) with time, keeping the characteristic form of the scaling spectrum. When 
the peak has disappeared, the cascade continues to lower x. causing a uniform, 
shape conserving, decrease of the occupations of the modes, and a flow of energy 
towards small x. 




(a) (b) 

Fig. 8. (Color online.) The function D(x,r) (Eq. (5.99)) at various times. Left panel: the 
filling of the modes, which proceeds till the disappearance of the leading particle peak. The values 
of t are, for the thick (blue) curves, from bottom to top: 0.1, 0.2, 0.3 (during this stage the 
leading particle acts as a source for soft gluon radiation), and for the thin (black) curves, from 
top to bottom: 0.5, 0.7, 0.9,1.0,1.1,1.2 (the leading parton has exhausted its energy and the peak 
has disappeared, while energy continues to flow to small x, the amount of energy in each mode 
decreasing exponentially fast). Right panel: energy is constantly injected into the system by a 
source located at x = 1 (see Eq. (5.106)). After a transitory regime, characterized by a uniform 
increase with time of the scaling spectrum, the system reaches a steady state. The values of r are, 
from bottom to top: 0.1,0.2,0.3, 0.4, 0.5, 0.6,0.7,0.8,0.9,1.0. 

We define 


T(x 0 ,t) = £( X q, t) = f dxD(x,r), (5.101) 

dT Jx 0 

where £(xq , t) is the amount of energy contained in the modes with x > Xq, and 
is the corresponding flux of energy, counted positively for energy moving 
to values of x smaller than Xq■ These quantities can be calculated explicitly. We 
have for instance 


£(x 0 i t) 


J dxD(x,r) = e 7TT erfc 


toq 

1 - x 0 


(5.102) 


with erfc(x) the complementary error function. We note that the fraction of the 
total energy “stored in the spectrum”, namely 


lim £( xq , t ) = e 7rr 


( 5 . 103 ) 
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decreases with time, and accordingly there is a non vanishing flux of energy reaching 
x = 0 


.F(0,t) = 2tt re'" . 


(5.104) 


It follows that the complete, energy conserving, solution involves a contribution 

D c (x) = n c (r)<5( x) with n c (r) = 1 — e~ nT , (5.105) 

somewhat analogous to a condensate where the radiated energy accumulates. Note 
that when r ~ 1/ y/w, corresponding to the disappearance of the leading particle 
into soft radiation, about 60% of the initial energy has flown into the condensate. 

It is interesting to consider also the situation where the leading particle is re¬ 
placed by a source that injects energy at a constant rate A at x = 1. In this case, 
we are led to look for the solution of the following equation 

d 

—D(x 7 t) = A5(l — x) +X[D], (5.106) 

where I[D] denotes the r.h.s. of Eq. (5.96). The exact solution of Eq. (5.106) with 
initial condition D(x , r = 0) = 0 reads 

D st (x,r) = A (5.107) 

27T Ak /a;(l — x) v ' 

and is plotted in Fig. 11, right panel. As time goes on, this solution converges to the 
stationary solution (A/2t:) / y/ x{l — x), keeping the shape of the small x spectrum, 
with just an overall time-dependent scaling. Remarkably, the small x scaling form 
D(x, t ) = f(t)/y/x is reached well before the stationary state is achieved. When the 
steady state is reached, all the energy provided by the source flows throughout the 
entire system towards the condensate at x = 0, that plays the role of a sink, while 
the population of the various modes stays unchanged. 


Let us contrast these properties with those of the DGLAP cascade. 95 A simpli¬ 
fied version of the corresponding evolution equation for the inclusive one particle 
distribution reads 

4 : D(x,t)=a[ d ~ . D(x/z,t) — a f -^D{x 7 t), (5.108) 

ot J x z(l - z) Jq 1 z 

where here the time variable is related to the virtuality Q 2 of the emitting parton, 
t = 1 uQ 2 /Qq. This equation has the form of Eq. (5.96) if one identifies 

—= a, and /C(z) = —r. (5.109) 

t*(x) 2(1 — z) 

Thus the DGLAP equation differs from the BDMPS-Z equation in two major as¬ 
pects. First, the kernel /C(z) is less singular near 2 = 0 and 2 = 1. Second, the 
rate of successive branchings is independent of the parent energy, i.e., it is constant 
along the cascade. We shall see that the latter property is what makes the major 
difference between the BDMPS-Z and the DGLAP cascades. 
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The solution of Eq. (5.108) is obtained using a Mellin transform .' 2 The initial 
condition reads D(v, 0) = 1, with D{v : 0) the Mellin transform of D{x , 0) = 5(1 — x). 
The solution can be expressed as the following integral 


D(x, t) 


pc-|-ioo 


Av 

2iri 


exp 


+ 7 ) t + pin 


1 

x 


(5.110) 


where (j>{y) stands for the digamma function and 7 ~ 0.58 stands for the Euler 
constant. 




Fig. 9. Model of a DGLAP cascade according to Eq. (5.108) (left panel), and with a source added 
at x = 1 (right panel). The curves corresponds to r = 0.2, 0.4, 0.6,0.8,1.0,1.2,1.4,1.6,1.8, 2.0 from 
bottom to top. In the right panel, the emergence of the scaling solution D(x) ~ 1/x is clearly 
visible, as well as the persistent deviation from it of the true solution at very small x. 


This integral representation allows us to verify a few properties that are relevant 
for our discussion. First, it is easy to check that the energy is conserved by the evo¬ 
lution: from Eq. (5.110) we have £(0,t) = f* AxD(x,t) = D(l,t) = 1. Furthermore, 
Eq. (5.110) has following asymptotic behavior when In 1/x t: 


D(x, t) 



(5.111) 


Here, in contrast to what happens with the BDMPS-Z cascade, the time depen¬ 
dence of the energy distribution does not factorize, and no obvious scaling behavior 
emerges at small x. The growth of the spectrum at small x is in fact tamed by 
the exponent in Eq. (5.111), making the spectrum integrable when In 1/x > t: all 
the energy remains in the spectrum, and no energy flows to x = 0 , unlike in the 
BDMPS-Z cascade. 

These features are illustrated in Fig. 9 where a numerical calculation of the 
DGLAP cascade is presented. One recognizes in the left panel features analogous 
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to that of the BDMPS-Z cascade, in particular the disaperence of the leading par¬ 
ticle peak, but no simple scaling pattern can be identified. In the right panel, the 
the solution in the presence of a source is displayed. The solution in that case 
approaches the stationary solution D(x) ~ A/x at not too small x, but it takes a 
long time to populate the very small x modes according to this 1/x law. In fact, for 
any time t , the true solution at very small x deviates from the stationary solution 
in the same fashion as in Eq. (5.111). 

The properties of the two cascades that we have discussed in this section can 
be understood in a broader context .' 2 The two cascades are governed by the same 
general equation, Eq. (5.96) with, in each case, a specific kernel K.(z) and a specific 
time scale t^(x). The kernel /C(z) controls how, in a given splitting, the energy 
is shared between the two offsprings. The time scale t*( x) controls the rate at 
which successive splittings occur. The characteristic features of the cascades result 
from two properties: energy conservation and the existence of approximate scaling 
solutions. The latter involves crucially t*(x). The role of the kernel KL{z) is less 
important. 

To better understand the role of t*( x), let us observe that Eq. (5.96) admits 
an approximate fixed point solution, which we refer to as a scaling solution, of the 
form 

D sc (x,t) = t —- (5.112) 

x 

Indeed, in the region where D(x,t) ~ D sc (x,t), there is complete cancellation 
between gain and loss terms. Note that this compensation occurs only for z > x 
in the integrals of Eq. (5.96), so that the scaling solution cannot be an exact fixed 
point of the equation. However the general solution is driven to this approximate 
fixed point. 

Now, since energy is conserved, and since it is continuously moving towards the 
lower values of x , the distribution will generically develop a divergent behavior at 
small x. Depending on the behavior of the scaling spectrum (5.112) when x —> 0, 
the approximate fixed point solution may lead to a divergent expression for the 
total energy £(0,f) contained in the modes between 0 and 1. In order for this not 
to happen, either the fixed point solution is never reached at very small x, or f*(;r) 
is an increasing function of x. The former occurs in the DGLAP cascade for which 
i*(x) = cst, while the latter is what occurs in the BDMPS-Z cascade for which 
f* ( x ) ~ yfx. Another way to see the role of the x dependence of f* (x) is to observe 
that in the DGLAP cascade, where f* is independent of x, it takes an infinite time 
to propagate a finite amount of energy from x = 1 to x = 0. In contrast, in the 
BDMPS-Z cascade, because the branching rate accelerates along the cascade (since 
i*(x) decreases with decreasing x (see Fig. 5 for an illustration)), the corresponding 
time is finite, and in fact of the order of t*. This is the reason why there is a finite 
flow of energy all the way down to x = 0 in the BDMPS-Z cascade while the energy 
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remains stored in the spectrum in the DGLAP cascade. For a more detail discussion 
see Ref.' 2 

At this point, it is useful to comment briefly on a peculiarity of the BDMPS-Z 
cascade, that makes the transition between the dilute, single-branching regime, and 
the the multiple branching regime completely smooth, with no sign of a change of 
regime as x crosses the value x s . Indeed, the perturbative solution, valid at r <C 1), 
is proportional to the BDMPS-Z spectrum 


D(x,t) 



(5.113) 


and this is already in the scaling form. We shall see that the scale x s becomes 
visible in the angular distribution, to which we now turn. 

The properties of the cascades depend also on the splitting kernel, that is, on the 
way the energy is distributed between the offsprings during a splitting. However, 
this turns out to have a minor effect on the main characteristics of the cascade, as 
compared to that of the transport time scale just mentioned, at least as long as 
the splitting kernel is not too singular. In that case the cascades develop as if the 
branching were completely democratic, with the two offsprings taking each half the 
energy of the parent gluon. The interactions responsible for the splittings can then 
be considered as local (in energy space). As already mentioned, these properties of 
the BDMPS-Z cascade that we have briefly listed, are typical of wave turbulence. 126 


5.2. Angular distribution 


We turn now the angular distribution. Our starting point is Eq. (5.90), that we 
write in Fourier space as 


(5.114) 


|D(x, u,t) =l~j d 2 JC(*) [f-D (* z , u ,t)-^=D (x, u, t) 


N c n 

4——cr(a :,u)D(x,u,t), 


-iO- 


with the Fourier transform given by 

/ d 

^D(x,0,f)< 

and similarly for a(x,u). We shall also consider the diffusion approximation ( 6' <C 
6) of Eq. (5.90) 


(5.115) 


^D(x,0,t) = J d zJC(z) 


D(x/z,0,t) D(x,6.t ) 

— z 


t*{x/z) 


t*(x) 


+ 


4 (xE) 2 \d0 


d 


— [qD{x,0,t)\ 


with q ~ 4:7ralN c n\n(6 2 /6 2 D ) and Op = mp /uj = mn/ixE). 


(5.116) 
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Solving Eq. (5.114) exactly is difficult. Part of the difficulty comes from the 
fact that the angular distribution has two distinct components: a hard component, 
corresponding to large angles produced by single hard scatterings, and a soft compo¬ 
nent that can be obtained as the solution of the diffusion equation (5.116). Each of 
these two components is strongly modified by gluon branching. The soft component 
admits moments, which is not the case for the hard component. The characteristic 
angle that marks the boundary between the soft and the hard components depends 
on x, i.e., on the amount of branching. It can be estimated by calculating the 
mean squared angle of the soft component. This provides a first orientation into 
the various regimes that characterize the general solution of Eq. (5.90), and that 
are illustrated in Fig. 10. 



Fig. 10. (Color online.) The various regimes of the physical processes (branching and scattering) 
that accompany the propagation of a fast parton with energy E >u c in a dense medium. The red 
curve stands for the root of (0 2 ). It separates the regime of single hard, large angle, scattering, 
from that of soft, small angle, multiple scattering. Accordingly, the soft component of the angular 
distribution lies below the red line, the hard component above. The line a 2 ui c ~ u s separates the 
region dominated by single branching from that of multiple branchings, u ~ E indicates the region 
where the leading particle propagates without splitting, and merely suffers momentum broadening 
due to its collisions with the constituents of the medium. 
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We define the typical squared angle as follows 

i r A^f) 

= (5 ’ U7) 

where D(x , 0, L) is here the solution of Eq. (5.116). This can be determined explic¬ 
itly in the three regimes characterizing the following ranges of x values: i) x < 1 , 
corresponding to the leading particle; ii) x s -C x <C 1 where single emission dom¬ 
inates; iii) x < x s , the regime of multiple emissions. We refer to Ref .' 0 for the 
details of the calculation, and discuss now the main results. 

When x < 1, 

(0 2 ) = § = 0 2 s (h L), d 2 (z, L) = g = ^. (5.118) 

This angle reflects directly the momentum broadening of the leading particle prop¬ 
agating in the medium over a distance L. 

In the regime x s -C x <C 1, where typically one emission occurs, the dominant 
contribution comes from momentum broadening of the radiated gluon from the time 
of its emission. The other contribution, corresponding to the momentum broadening 
of the leading particle before the splitting, is suppressed by a factor x 2 . Keeping 
the dominant contribution, we get 

{6 2 )= l -0 2 s {x,L). (5.119) 

The factor 1/2 originates from the average over the emission time. 

Note that for realistic situations co c < E, the maximal frequency is thus x c = 
uj c /E <1. In this case, x s <C x <C x c , and hence, these rare radiations are on 
average confined to the angular corona, 

0 C « « « -4 d c . (5.120) 

a i 

where the critical angle 

e c = (gi 3 )- 1/4 , (5.121) 

corresponds to the smallest BDMPS-Z radiation angle. The third regime, x x s , 
is dominated by multiple branchings at parametrically large angles, 

0 > \o c . (5.122) 

It requires a non perturbative treatment and the result is (see Refs . 6 '’ 70 for details 
and Ref . 68 for a similar discussion) 

io^\ = _l r g 

' ; 4a l(xE) 3 

As was the case for the previous regime, this angle can be identified with that 
corresponding to the momentum broadening of the observed gluon since the time of 


= 0l(x). 


(5.123) 
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its emission. The major difference with the previous regime is that now 9*(x) does 
not depend on the size of the medium. This is because, in the multiple branching 
regime, only gluons that are emitted off the leading particle at a distance, that is 
shorter than their characteristic stopping time from the medium boundary 

L, escape the medium. Accordingly, their measured squared angle is typically 
0 2 {x) ~ qt*(x)/(xE ) 2 , which recalling Eq. (5.94) agrees with Eq. (5.123). 

The complete distribution in the three regimes can also be determined explicitly. 
Again, we refer to Ref. 70 for details. 

When x < 1, the distribution takes the factorized form 


D{x,O,L)~V(0,L)D(x,L), 


(5.124) 


where D(x,L) is the energy distribution given by Eq. (5.99) (for x < 1). In the 
regime of multiple scatterings, i.e., for 9 2 < 9 2 (L,x) 


V(0,L) 


47T 


exp 


e 2 ' 

. 0 2 (x,L)_ ' 


(5.125) 


In the opposite case, 9 9 s {x,L), a single hard scatterings dominates the distri¬ 

bution and we have 


V(9,L) 


nN c g A L 

(. xE ) 2 6> 4 ' 


(5.126) 


In the second regime the energy distribution is given by the leading order 
BDMPS-Z distribution, 

D{x, L) ~ —, (5.127) 

and we expect also the angular distribution to be given by a single radiation that 
undergoes multiple scatterings. For sufficiently small x, as already emphasized, the 
angular deviation of the radiated gluon is larger than that of the leading parton. 
Accordingly one can neglect the momentum broadening before the emission. We 
get, 37 

r L dt 

D(x,0,L)~ -V(x,0,L-t)xlC(x). (5.128) 

Jo t* 

Finally, in the fully non-perturbative regime, i.e., x x s , and in the soft region, 
both multiple scatterings and multiple branchings must be resummed. A solution 
can be obtained as a power series in the number of scatterings.' 0 This is achieved 
by using an approximation that is equivalent to the diffusion approximation, hence 
not accurate for the very the tail of the distribution (see below), but it allows us to 
determine the distortion of the main peak. The result reads 


4tt / 0 2 \ 


D(x,0,L) 


D(x,L ), 


(5.129) 
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where the entire dependence on 0 is in the scaling function 77 , given by 


77 ( 2 ) 



da J 0 (2y/za) ^ c n (-a 2 )" 

n—0 


(5.130) 


where Jo is a Bessel function, and we have the normalization property. 

dz r)(z ) = 1. (5.131) 

In Fig. 11 we have plotted the angular distribution 77 ( 2 ) in the multiple branching 
regime. For the numerical evaluation we have computed the first 500 terms in the 
multiple scattering series (5.130). 70 

Consider finally the single scattering limit achieved when 6 6^{x). This is 

easily found to be 

D (x, 0. L) « N ^y ( ^ D{x, L), (5.132) 

where N c ng 4 t*(x) / (2xE) 2 ~ 0 2 (;r) is the typical angle squared of gluons in the 
multiple branching regime. 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 


V~z 


Fig. 11. (Color online.) The angular distribution (5.130) of a gluon in the multiple branching 
regime (solid black line) compared to an exponential distribution with identical first two moments 
(dashed purple line). Here z = 0' 2 jO'l(x). 


5.3. Application: Energy imbalance in dijet events 

As we have argued in this section, a striking feature of the BDMPS cascade is that 
it provides an efficient mechanism to transport energy down to the lowest accessible 
gluon frequencies. Since soft gluons are typically emitted at large angles, the cascade 
naturally populates soft gluon modes at large angles with respect to the jet axis. 
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We end this section with a brief analysis of the phenomenological implications of 
this mechanism. 

Recently, the CMS collaborations investigated the energy imbalance in strongly 
asymmetric dijet events. 18,26 As they are born back-to-back jets have roughly equal 
energy and momenta. The observed asymmetry is due to the fact that the leading 
jet (the most energetic of the two) traverses a shorter distance in the medium than 
the subleading one, and hence loses less energy. Surprisingly, the energy balance is 
recovered at quite large angles, R > 0.8, and most of the “missing” energy is carried 
by very soft particles. 18 11 



Fig. 12. The difference of the angular and energy distributions of sub-leading and leading jets for 
ojbh = 0.5 GeV. The histograms account for the four binnings of energies: [0-1] GeV (grey), [1-2] 
GeV (yellow), [2-4] GeV(orange), [4-8] GeV(green), [8-100] GeV(red). In the first angular bin we 
observe a large imbalance of energy in the hard particles. This energy is partly recovered at large 
angles by very soft particles. The cumulative energy is given by the full line. 

As a very rough model for an asymmetric dijet event, we consider a leading jet 
and a subheading one that have traversed, respectively, the path lengths L = 5 
fm and L = 1 fm in opposite directions. We use the following set of parameters: 
q = 1 GeV 2 /fm, a = 0.3, and E = E 1 = E 2 = 100 GeV. First, we find that the 
broadening of the subleading jet, 0 S ~ 0.02 (cf. Eq. (5.118)) is negligible compared 
to the opening angle of jet 0j e t = 0.3, as observed in the data. 

We evaluate the energy imbalance by integrating the difference D 2 (x,6) — 
D i(cc, 0), where D\ and D 2 are the medium-induced gluon distributions of the lead¬ 
ing and sub-leading jet respectively, in the different bins in energy and angle. The 
minus sign accounts for the fact that the constituents of the subleading jet prop¬ 
agate in opposite direction to those of the leading jet. The angular distribution 
is assumed to be Gaussian with a width determined by an interpolation between 


n To our knowledge, this observable was only addressed in the context of the Q-PYTHIA Event 
Generator 127 which reproduces the qualitative features seen in the data. 
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the mean squared angles in the three regimes discussed above. Finally, the evolu¬ 
tion equation Eq. (5.96) is solved numerically where soft emissions with frequencies 
u) < wbh are cut-off to regulate the unphysical soft singularity in the branching 
kernel (3.51). 67 

In Fig. 5.3, is plotted, in the form of histograms, the energy imbalance integrated 
on various angular and energy intervals. 6 ' At small angle 9 < 0.2 there is a deficit of 
hard particles with energies ui > 8 GeV (red negative histogram). As we increase the 
angle the energy imbalance turns in favor of the subleading jet and it is dominated 
by very soft particles with energies ui < 1 GeV (positive grey histograms). The 
full line corresponds to the integrated energy imbalance for angles from 0 to 6. 
The imbalance is expected to be recovered at large 0, when 9 —> n/2. The energy 
transported at large angles is clearly dominated by very soft gluons °. The jet 
structure that emerges from this plot is in semi-quantitative agreement with the 
recently released data by CMS. 26 

6. Decoherence of the jet core 

In this section we discuss briefly how the two cascades, vacuum and medium in¬ 
duced, can interfere in a jet produced in heavy ion collisions. The issue we want 
to address is to what extent color coherence, that constrains the pattern of succes¬ 
sive parton branchings in vacuum, is altered when the jet propagates in a medium. 
As a paradigmatic example, we consider the radiation (in vacuum or in a QCD 
medium) of a quark-antiquark pair in the singlet state (the generalization to dif¬ 
ferent color configurations is straightforward). This pair, which we refer to as an 
antenna, can be thought of as being produced by the decay of the highly virtual 
primary parton. 59 ~ 64 

In vacuum, the radiated spectrum off the quark reads 3 (see discussion in Section 


2 . 1 ) 



(6.133) 


and similarly for the antiquark, i.e., dP = d P gq + d P gq . Here, the theta function 
accounts on average for the suppression of large angle radiation (9 9 qq ) caused by 
color coherence. The total radiation spectrum exhibits two radiation cones aligned 
on the quark and the antiquark momenta, respectively. Eq. (6.133) is at the basis 
of Eq. (2.2). 

The probability for the quark of the antenna to radiate soft gluons after travel's- 


°These features are absent in the angular distribution of primary gluon radiations, as shown in 
the early work, 128 confirming the important of multiple branchings to transport energy at large 
angles. 
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ing a QCD medium is given by, 59,60 
ol s Cf 


d Pgq — 


0(9 q q - 0) + A m ed 0(0 - dqq) 


■ dd 2 do; 


(6.134) 


where 

r X 

ill 


A 


med — 


l 

N? - 1 


(Tr Up, 


1 — exp 


(6.135) 


and U is a Wilson line in the adjoint representation, evaluated along the trajectory of 
the quark, and the harmonic approximation has been used to evaluate the correlator 
of the Wilson lines. Here rj_ = 9 qq L stands for the antenna transverse size and 


l± = 



(6.136) 


is a length characterizing the range of coherence in the transverse plane. The 
quantity A me( j may indeed be interpreted as a decoherence parameter: when rj_ <C 
Zj_ , or equivalently 0 C Oj et , where the BDMPS-Z critical angle 9 C is defined in 
Eq. (5.121), we have A med ~ r\H\ -C 1 and we recover Eq. (6.133). In this regime 
the inner structure of the antenna is not resolved by the medium and the antenna 
radiates coherently as in vacuum. 



Fig. 13. The soft gluon emission spectrum off the quark constituent of a singlet antenna with 
opening angle 6 qq = 0.2, according to Eq. (6.134), in the presence of a medium with A me d = 0.5 
(solid line). Here = asCp/n = 1. On average vacuum radiation is confined within 6 < 6 q q, while 
the medium-induced radiation is radiated at 6 > 6 qq . The limit of opaque medium, given by 
A me d = 1, is marked by the dashed line. Figure taken from Ref. 60 


On the other hand, the case rj_ l± ( 6 C <C Oj e t), corresponds to full decoher¬ 
ence of the antenna, and we have A med — 1. Here, the angular constraint disappears, 
















March 23, 2015 0:17 


World Scientific Review Volume - 9.75in x 6.5in 


QGP5-BM-FINAL 


46 


J. -P. Blaizot and Y. Mehtar-Tani 


and the quark and the antiquark radiate independently from one another 60p , 


d P, 


91 


a s C F d 8 2 dw 
7T 9 2 U) 


(6.137) 


This discussion applies also to gluons radiated inside the medium . 61 ~ 64 The deco¬ 
herence parameter given by (6.135) can be rewritten as A me <j ~ 1 — exp(L/idecoh), 
from which we extract the decoherence time, 


blecoh — 


1 

qOl 


ii 


1/3 


(6.138) 


The antenna spectrum (6.134) is plotted in Fig. (13) for A me d = 0.5. 


The antenna setup serves as guideline to discuss the decoherence that a jet 
experiences when passing through a dense medium. Consider a jet with opening 
angle Bj e t- While inclusive jet observables in vacuum are characterized by two 
scales, Q = 0j e t^, that correspond to the jet transverse mass typically of the order 
of 100 GeV, and Q o ~ Aqcd, in the presence of a dense medium, an additional 
(hard) scale comes into play, Q me d = max (rj 1 , y/qL) , where here rj_ = 0j et L. 
This medium-generated scale is of the order of a few GeV’s, therefore, we typically 
have Qmed Q■ The medium affects the collimated parton shower, and thus 
causes its decoherence, only at scales of the order or smaller than Q me d- Hence, the 
medium alters the property of angular ordering by allowing for additional soft gluon 
radiations up to the limiting angle Qmed / W - 63 ’ 64 As a result of this aforementioned 
separation of scales a large portion of the jet is barely affected by its interactions 
with the medium. 





Fig. 14. An illustration of a quark-antiquark antenna of typical transverse size r± = Q qq L in a 
medium of length L, in the coherent regime (left) where the typical in-medium resolution length 
l± = ( qL)~ 1 / 2 r±. In this case, the medium doest not resolve the inner structure of the 

antenna, and thus interacts with its total charge. In the opposite case (right panel), l± =<C r±, 
the medium resolves the quark and antiquark color charges, causing their decoherence. Then each 
of the antenna constituents radiates independently. 


Comparing the two medium scales one is left with two regimes: 


P A analogous phenomena occurs between initial and final state radiation in a dense medium . 129,130 
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(1) r_L -C l_ l (or 0j e t <C 9 C ): in this situation the jet is on average not resolved by 
the medium and remains color coherent. Multiple interactions with the medium 
will induced gluon cascades off the total charge, typically at large angles (cf. 
Eq. (5.122)), 

jot- (6.139) 

°^S 

The unresolved vacuum cascade looks as if it developed outside the medium. 

(2) In the opposite case, r± >> l± (or 0j et 9 C ), the medium resolves more charges 
inside the jet, which results in an enhancement of the medium-induced radiation 
rate. In other words, the medium resolves a certain number of sub-structures 
characterized by the opening angle 9 C , that evolve independently of each other. 
On the other hand, the inner structure of these sub-jets remain unresolved by 
the medium. 131 Consequently, the unresolved system loses energy coherently 
with a rate proportional to the total charge contained within 9 C . 

For realistic situations it was noted 131 that a core containing most of the jet 
energy remains unresolved by the medium. Secondary substructures carry typically 
smaller energy fractions and are expected to be influenced significantly by medium 
effects. Under these circumstances, an unresolved jet lose energy as a single parton 
whose final distribution is then given by the rate equation (5.96). In this approxi¬ 
mation, we can write the jet spectrum in nucleus-nucleus collision as a convolution 
of the hard parton spectrum due to incoherent binary collisions, and the probability 
distribution to find a parton with energy p± that originates from the nascent parton 
with momentum p' ± = p±/x, 

dM A /' 1 dx dJVS 

—U— = N coll — D(x,p±/x)-^— (p±/x) (6.140) 

d 2 p± J o x d^pj_ 

where the distribution D(x,p' ± ), was discussed in Section 5.1, identifying p' ± = E 
and stands for the number of binary collisions. It is obtained by solving nu¬ 
merically Eq. (5.96), where finite size effects are taken into account (cf. Eq. (2.6)) 

. This procedure is analogous to that used to computed quenching of hadron spec¬ 
tra, 40,42 ’ 132 but proves to be best justified for high-p* jet spectra. The resulting 
nuclear modification factor, Raa = -^randiVj^/dfVjgj, is plotted in Fig. 15(a), 133 
where the uncertainty bands shows the sensitivity of this calculation to variations 
of the medium length L and the in-medium diffusion coefficient q, via the variable 
w c = 60 — 100 GeV. In addition the sensitivity to the soft scale wbh = 0.5 — 1.5 
GeV is also shown, is implemented as a cut-off of the unphysical singularity in the 
branching kernel (3.51). For q = 2 GeV 2 /fm and L = 4 fm, we get oj c = 80 GeV 
which corresponds to the average value. These values are consistent with the usual 
expectations. 123 

In the CMS analysis the jet opening angle is chosen to be 0j et = 0.3, it follows 
that rj_ ~ 1.5 fm and l± ~ 1.1 fm. Obviously, the jet is neither in the fully coher¬ 
ent regime nor in the decoherent one. Hence, decolrerent radiation off secondary 
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substructures may be important. It was suggested in Ref. 133 that a signature of 
partial decoherence may be found in intrajet structures, e.g., in the fragmentation 
functions, or equivalently the energy distribution of particles inside the jet-cone, 
D(x) = xdN/dx. 


05 



50 100 150 200 250 300 350 400 

P± [GeV] 



l 


(a) 


(b) 


Fig. 15. (a) Calculation of the nuclear modification factor with oj c = 80 GeV as a function of 

jet p± for central Pb-Pb collisions. The solid (red) band includes the variation of ubh around a 
central value of 1.5 GeV. The dashed (grey) band includes, in addition, a variation of u c S [60, 100] 
GeV. The experimental data are taken from." 4 (b) Upper panel: the longitudinal fragmentation 
function plotted as a function of £ = Inl/x. Lower panel: the ratio of medium-modified and 
vacuum fragmentation functions. The experimental data are taken from. 134 See text for further 
details. 


The strategy to identify such effects consists in looking for departures from the 
fully coherent limit. In the coherence regime, the unresolved vacuum shower loses 
energy out of the jet-cone as a single parton. Hence, we can effectively factor out 
the in-medium distribution D mex i(z.p±) described by Eq. (5.96) from the vacuum 
evolution (see Ref. 133 for details): 

r 1 d~ 

D coh (x)= / —D vac (x/z,Q)D med (z, P± ), (6.141) 

where D vac (x/z,Q) is solution to the MLLA evolution equation (2.2) with Q = 
0jetP_L- The result is plotted in Fig. 15(b) and is given by the grey band. Eq. (6.141) 
properly matches the pp baseline with D me d{z,p±) = d(l — z) (see the blue curve 
in Fig. 15(b)). The semi hard sector is well described and is consistent with the 
nuclear modification factor above. However, there is an excess of few soft particles 
that is not accounted for. We include now the decoherent contribution, AD decoh 
which is computed to double logarithmic accuracy by convoluting the probability 
to form gluon-quark antenna in the medium which then subsequently radiates de¬ 
coherent soft gluons according to Eq. (6.134) (see Ref. 133 for details). The total 
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distribution, D tot ~ D coh + A D decoh is displayed as red band in Fig. 6.134. This 
in-cone decoherent radiation proves to be sufficient to account for the excess of soft 
gluons seen in the data. 


7. Conclusions and outlook 

In this review, we have analyzed the dominant features of the medium-induced 
parton cascade, relying on the BDMPS-Z mechanism for medium induced gluon ra¬ 
diation. The medium induced cascade exhibits remarkable properties that originate 
on the one hand from the specific mechanism of medium induced radiation, on the 
other hand from the rapid loss of coherence among the radiated gluons caused by 
their multiple interactions with the medium constituents. These properties manifest 
themselves in a specific mechanism for the transport of energy along the cascade, 
with a characteristic angular structure. The transport of energy involves a constant 
flow down to the lowest accessible frequencies. The angular structure is such that 
the softest quanta are emitted at large angles. We emphasize that these features, 
reminiscent of those revealed by the analysis of recent CMS data on the angular 
distribution of the missing energy in inbalanced clijet events, are intrinsic proper¬ 
ties of the BDMPS-Z cascade. They do not result from a specific coupling to the 
surrounding medium, nor are they influenced by the collective dynamics of this 
medium. In fact the only dependence on the medium is hidden in the value of the 
jet quenching parameter q. 

Much remains to be understood about the BDMPS cascade. In particular the 
treatment of what happens at the lower end of the cascade, where gluon frequencies 
become smaller than ubh needs to be improved. In view also of the potential role of 
this cascade in the thermalization of the quark-gluon plasma, 125,139 understanding 
how energy carried by the soft radiation is dissipated in the medium is an important 
task. An important ingredient entering the description of the medium induced 
cascade is the jet quenching parameters, q. We have seen that this coefficient 
receives sizable radiative corrections. The calculations of these corrections is at 
present limited to the most singular contributions (double logarithmic). We have 
argued that it is the singular nature of the correction that allows us to absorb it 
in a redefinition of q. It is unclear whether the simple description of energy loss 
and momentum broadening as local transport phenomena would hold in a more 
accurate treatment of these radiative corrections. 

Aside from the medium-induced cascade that we have just discussed, a jet in 
a medium is accompanied by another, vacuum-like, cascade. The two cascades 
are geometrically separated. The vacuum-like cascade is a collimated and coherent 
cascade, constituting the jet core, characterized by collinear and angular ordered 
splittings triggered by the initial hard collision. The medium-induced cascade devel¬ 
ops at parametrically large angles by successive branchings of the primary radiated 
gluons off the unresolved total color charge of the jet core. Much remains to be 
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understood about the interaction between the vacuum cascade and the vacuum in¬ 
duced cascade. An interesting phenomenon takes place when the medium charges 
resolve substructures in the inner jet core. The coherence of the resolved charges 
is broken, allowing them to radiate independently in the soft sector. This mecha¬ 
nism for medium-induced soft radiation is of a different nature than the BDMPS-Z 
mechanism. It causes the decoherence of vacuum radiation which is not broadened 
by the medium and hence is expected to remain collimated. This mechanism may 
provide a natural explanation for the excess of soft particles measured in the jet 
fragmentation region. Finally, in addition to the effects on the intrajet structure, 
color decoherence is expected to be significant for processes at NLO, e.g., observ¬ 
ables involving multi-jet events. 

A complete description of parton cascades in a QCD medium goes along with 
a refined understanding and modeling of the medium dynamics. 92,135-13 ' This, 
together with the details of the experimental setup, requires the development of 
event generators, 92,138 an effort undertaken by many groups. 85-89,92 Clearly, the 
microscopic physics that has been discussed in this review needs to be implemented 
in these generators. While this is the case, to some extent, of the BDMPS-Z cas¬ 
cade, 85,88,89 as far as we know, the decoherence mechanism as well as the renormal¬ 
ization of the quenching parameter are not yet accounted for. 
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